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1. Introduction. We consider those properties of turbulence in fluid of zero mean 
motion which can be deduced from the assumptions of spatial homogeneity and isotropy. 
The 7-components of the velocity fluctuations at two points P(x) and P’(x’) will be 
written as u; and u’ . The spatial separation of P and P’ is denoted by the vector 


of magnitude r, where 


and repeated indices imply summation over the values 1, 2 and 3. 

It may be shown without difficulty’’*® that the condition of isotropy requires the 
correlation between u,; and u} , where 7 and j have arbitrary values, to depend only on 
the geometrical configuration defined by — and unit vectors in the 7- and j-directions, 
and a single scalar function of r*. It can be represented as the typical component of a 
tensor of the second rank in the following manner: 


Pe EP on el _1 of 1, 2f) | 
uwus = Ri = u | Or ar esi + (; + 2 ar b;; (1.1) 





where the over-bar denotes a spatial mean value, f(r) is the scalar function, and 


5,;; is the unit tensor whose value is 1 if 7 = j, and 0 otherwise. The triple correlation 
between velocity components-at P and P’ has a similar dependence on a single scalar 


function; 


uuu, = T%; 
_ | (k= 2 ok /ar t+ Hy ak/or) _k | 
= u ( 2, Jest + ( Or (Gj: + Ses€;) Or b;€% |. (1.2) 








*Received May 3, 1947. 
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The scalar function k(r) is in this case an odd function of r, of order 7° when r is small. 
It is clear from (1.1) and (1.2) that f(r) and k(r) are the correlation coefficients for 
particular values of the indices and the é; , viz. 


u’-f(r) = (us); , 


u!®-k(r) 


(wu jury) ke 


eter 


All higher order correlations which involve at least one velocity component at each of 
the points P and P’ depend on more than one scalar function of r. 

The equations of motion of a viscous incompressible fluid contain both linear and 
quadratic terms in the velocity components, and it is consequently possible to relate 
the double and triple velocity correlations. Expressed in terms of the scalar functions 
f(r) and k(r), this relation becomes the equation for the propagation of the double 
velocity correlation and has the form 


— 
 - u(e + i.) + ayu" 7” + 4 9), (1.3) 
where dashes to f and k denote differentiation with respect to r, and ¢ is the time. Equa- 
tion (1.3) has a simple and useful form for the particular value r = 0, when it describes 
the rate of decay of energy of the turbulence; 
72 72 

an ~ ee «+ (1.4) 
where } is the length parameter previously introduced by Taylor.* 

Since the difficulties of measurement become very great as the order of the correla- 
tion increases, it is inevitable that Eq. (1.3) should occupy an important place in any 
practical theories of turbulence. The purpose of this paper is to extract as much in- 
formation from it as is possible with a minimum of further assumptions, and in particular 
to deduce the rates of energy decay which are consistent with certain types of behaviour 
of the function f(r). It is important to appreciate that Eq. (1.3) represents all the in- 
formation about the function f(r) which we possess, and is clearly insufficient to permit 
a complete solution to be obtained. This lack of information for the determination of 
correlation functions is inherent in averaged equations, and is the penalty paid for this 
use of the statistical method. Our plan is therefore (a) to discuss suitable limiting cases 
for which a solution of (1.3) is possible, and (b) to discuss the consequences of simple 
hypotheses so that their validity can be put to the test of experiment. Types of solution 
for which the function f is “self-preserving”’, i.e. for which f is a function only of r/L, 
where L is a length which depends on the time ¢, were introduced by v. Karman and 
Howarth.’ The consequent gain in simplicity of the mathematics is considerable and 
it is with such solutions that we shall largely be concerned here. There is available 
sufficient experimental evidence to indicate that the theoretical solutions: correspond 
in some respects with reality and are worth pursuing. 

2. Loitsiansky’s invariant. We shall have need later of a simple deduction from the 


3G. I. Taylor, Statistical theory of turbulence, Proc. Roy. Soc. (A) 151, 421-478 (1935). 
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basic Eq. (1.3), which was pointed out first by Loitsiansky.* Multiplying both sides of 
the equation by r* and integrating over the range 0 to ~, we find 


8 (ie. [os ) = 
2 (u [orf ar = 0 


provided that r*f’ — 0 and r*k +0 as r— o. Physically it seems reasonable to suppose 


that the convergence conditions are satisfied; then 


ul. | r*f dr = A, (2.1) 


“0 


where A is a constant during the decay process. 

Loitsiansky remarks that the relation (2.1) and Eq. (1.3) are consistent with an 
analogy between the propagation of fu’’, and the propagation of heat in a spherically- 
symmetrical five-dimensional field. In such a space the Laplacian operator has the form 


so that the last term of Eq. (1.3) represents the effect of molecular conduction provided 
that fu’* is the analogue of temperature. The first term on the right side of (1.3) must 
be interpreted as the effect of convection, and (2.1) shows that this convective effect 
is such as to leave the total quantity of heat constant. On the basis of this analogy 
Loitsiansky describes the constant A as a measure of the total quantity of disturbance 
to the fluid, which is uniquely determined by the initial conditions of the turbulence. 

It is not without interest to notice the behaviour of the general moment of the 


function f(r). Thus 


a! wa 13 
ry (u | rf ar) = (4 — m)u | 


poo 


kr"~* dr + 2(m — 1)(m — 4)vu”” l r"~?f dr (2.2) 


“0 


provided that m > 1 and (r”"'f),. = (r"k)~/ = 0. When m = 1 or 0, special formulae 


are necessary, Viz. 


S (w?- f rf ar) = au. | k dr — 6vu’’ , (2.3) 
< (w°. | f ar) = 4u’- | * ar + 8vu”” | _ dr. (2.4) 


Now all the experimental evidence suggests that the function f(r) is everywhere positive 
and monotonic decreasing, and the signs of the integrals containing f and f’ can be 
predicted with safety. There is not very much data about k(r) but measurements at 
the Cavendish Laboratory show it to be everywhere negative (or zero). It must certainly 
be negative for small values of r in order to give a positive contribution to the rate of 
change of mean square vorticity due to random extension of the vortex lines.* If the 





‘L. G. Loitsiansky, Some basic laws of isotropic turbulent flow, Central Aero- and Hydro-dynamic 
Institute, Moscow, Report No. 440, 1939; also N. A. C. A. Tech. Memo. 1079. 
6G. K. Batchelor and A. A. Townsend, Decay of vorticity in isotropic turbulence, Proc. Roy. Soc. (A) 


190, 534-550 (1947). 
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signs of the integrals involving k are predicted on this basis, Eqs. (2.2), (2.3) and (2.4) 


show that 
Ld (w. [ry ar) se. woe 


“0 


a> 2 weed (2.5) 





<0 if m = 3, 2,1, 0. 


These results have a bearing on the behaviour, during decay, of the spectrum func- 
tion describing the spatial structure of the energy of the turbulence field. For if u’?- F(u) du 
is the amount of energy lying within a small range about the wave number uz the functions 
f(r) and F(u)/2(2r)’” are Fourier transforms,° and 


F(u) = 4 | cos 2rpr- f(r) dr. (2.6) 
0 
The relations (2.5) thus become, when m is even, 
a 2F) > ’ 
ry | w { _ = 0 if m 


If the energy spectrum function is expanded in powers of yu’, viz 


4!” ea 72 du F) eB (ae) ams 
‘F(u) u FO +5 x (a iT 4l au! Tins ’ 


(2.7) 


AllV 
ce 








then the effect of decay is to decrease the coefficients of »’° and y’, to leave the coefficient _ 
of u* constant and to increase all other coefficients. 

3. Self-preserving solutions at Reynolds number which are not large. By Reynolds 
number which are not large is meant a state of affairs in which d is not small compared 
with other lengths associated with the function f(r). Since the length \ is already present 
in the basic Eq. (1.3) (in view of the expression for du’*/dt), it will be convenient mathe- 
matically to use \ as the scale factor L in any solution for f(r) which preserves its shape 
over a range of r which includes small values of r. Several possibilities can be considered, 
the simplest of which is that the correlation functions preserve their shape for all values 
of r. 

The hypothesis whose consequences are to be examined is that 


fr =f), kr) = kW (3.1) 


for all values of r, where y = r/X. Equation (1.3) then becomes 
67 AD yt 1va (me , M8), (EL 4 4 as) 
of — Ee te 2» \dy’ yv/* \ayt vay 
or, in terms of the number R, = u’d/v, 


af ,4df , 5ydf 5f) ? dk, (, af (2 4 - 
(a + yayt 2dp* !) t+ op, at (v Ht) +5 + 5% \ay t 7] =o C4) 








*G. I. Taylor, The spectrum of turbulence, Proc. Roy. Soc. (A) 164, 476-490 (1937). 
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Each of the expressions within circular brackets is a function of y only, while the co- 


efficients outside the brackets depend only on t. 
It is worth noting that the second coefficient, viz. \” dR,/2vR, dt, is a constant for 
a very general class of energy decay law. If u’® decays as some power of ¢, the energy 


equation (1.4) shows that the decay laws will be 


w?n~t, N=  , Ry~ t™ (3.3) 
and hence 
? dR _ 5(1—n) | 
QR, dt 2 (3.4) 


An exponential decay of u’’ also makes this factor constant. When the second coefficient 
is constant, the first two groups of terms in (3.2) are functions of y only and the equation 
can only be satisfied by 

R, = constant. 


This leads to n = 1 provided the constant value of R, is not zero (in which case we should 
have n > 1, i = @). When the energy decay law is not such as to make )d” dR,/2vR, dt 
constant, there is another way in which Eq. (3.2) can be satisfied. This alternative 
method is suggested by the work of Sedov,’ and will be discussed in section 5. 

Now the law of energy decay is already fully determined by the assumption of self- 
preservation of the correlation function f for all values of r. For in this case the condition 


(2.1) becomes ‘ 
wr [VS av = 


provided that the integral converges, and thence u’*\° is constant during the decay. 
The energy equation then gives the decay laws as 
-2 5/2 + - 
yl ae AT. ”? = Apt, Rx = A _ (3.5) 


where A is a constant and ¢ is measured from the instant at which 1/u’ = \ = 0. This 
is a power law of energy decay (with n = 5/2) and as shown above, Eq. (3.2) can only 
be satisfied when R, is constant. Since R, also varies with ¢ according to (3.5), the two 
requirements can only be consistent if 


’ 


Under these circumstances Eq. (3.2) becomes 


dj , af (4 ) ‘ 
dy’ dy + v + 5f bt 0 » (3.6) 
and the solution which makes f = 1 whenr = 0 is 
fw =e". (3.7) 


7L. I. Sedov, Decay of isotropic turbulent motions of an incompressible fluid, C. R. Acad. Sci. U. R. 8. 8. 
42, 116-119 (1944). 
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The requirement R, = 0 shows that the triple correlations have no influence on the 
double correlations under the conditions for which a completely self-preserving solution 
is possible. 

Thus it has been shown that a solution in which the function f is completely self- 
preserving is only possible at large decay times and is described by (3.5) and (3.7). 
This suggests that we should examine conditions as t > @ in order to see if such a self- 
preserving solution does in fact exist there. It is not difficult to see at once that the 
answer is likely to be in the affirmative. The correlation function f(y) certainly has 
the same (parabolic) form for all values of ¢ for y < 1, and the only alternative to an 
approach to a definite shape as t > © is an oscillation of the remaining part of the 
curve. Such an oscillation does not seem appropriate to the problem. 

However, more definite evidence that f(y) is independent of t when ¢ is large can be 
obtained from the basic equation. For consider (2.3) in the form 


oa [fs r) - Sup Pe: ee 
2 (Re Pi fas)= mf] kat —6e. 





All the measurements of f(r) which have hitherto been made have shown it to be every- 
where positive and we may safely assume the expression within brackets on the left 
side to be positive and, in view of the existence of a parabolic variation for r/A < 1, 
non-zero, for all values of ¢. Also, as discussed in section 2, the evidence is that k is 
everywhere negative (or zero), so that 
po 72 

g(a 5S ax) < 0% as 
Suppose now that when ¢ is large, the decay of u’ and \ conforms to the general laws 
(3.3). In the first place, n cannot be less than or equal to one, for (3.8) then requires 
Rx fo r/X f d r/d to decrease indefinitely and to become negative as t > ~. Secondly, 
the inequality (3.8) requires the expression in brackets on the left side to vary as some 
power of ¢ not greater than 1 — n, whereas the factors in this expression show that it 
varies as some power of ¢ not less than 1 — n. Hence the power can only be 1 — n and 
Jo r/A f d r/d tends to a constant when ¢ is large, which shows that f tends to become 
a function only of r/X. Notice that the assumption k < 0 is over-sufficient; the deduction 
is only rendered invalid if fy k d r/d is positive and increases as some power of ¢ not less 
than (n — 1)/2. 

It is thus possible to make quite definite predictions about the turbulence when ¢ 
is large. Making assumptions about f and k which are well supported by experiment, 
and assuming that the energy decay follows a power law (3.3), it can be shown that 
f(r) preserves its shape when ¢ is large. Then Loitsiansky’s invariant relation shows that 
the decay laws must be as in (3.5). Finally the fundamental equation for f gives the 
solution (3.7). The.dependence of these deductions on a decay law of the type specified 
by (3.3) when ¢ is large does not seem likely to be critical. 

4. Solutions obtained by neglecting the triple correlation. Two Russian authors, 
Loitsiansky* and Millionshtchikov,*® have each discussed the solutions of Eq. (1.3) 
which are obtained when the term describing the effect of the triple velocity correlation 


8M. Millionshtchikov, Decay of homogeneous isotropic turbulence in a viscous incompressible fluid, 
C. R. Acad. Sci. U. R. 8. S. 22, 231-237 (1939). ° 
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is ignored. Their work is an extension of the “small Reynolds number’”’ solution first put 
forward by v. Kaérman and Howarth.’ There is considerable indirect evidence for the 
belief that neglect of the triple correlations is only permissible at a late stage in the 
decay of turbulence so that the resulting solutions ought to be compared with that 
deduced in the previous section for ¢ large. 
The equation to be solved is 
eh = 2, (Put 4 Set), (4.1) 
ot or r or 

which is similar to the equation for the propagation of heat in a spherically-symmetrical 
five-dimensional field in which there is no convection. Using this analogy the solution 


is known to be 
u"f = ; — fff F(s, O)e~?’’*"' da, dx. dx, dx, dx; (4.2) 
‘ (Savt)’’* Je ed , 


where 
F(r, ) = w’’f, 
‘ ; 5 5 
P= LeE—-z), P= LE, #®= Dae. 
n=1 1 1 
The law of energy decay is obtained by putting r = 0 in (4.2), i.e. 


ae F(s, 0) e*'’*"' da, da, dx, dx, dx 
— il ) + dos de 





, ' an 
= ear ge F ; 0 — a Is. 4.3 
piers, lll - 


Thus, in general, the solution depends upon, and is uniquely determined by, the function 
F(r, 0). 
If the initial state of the turbulence is such that 


F(r, 0) = 0, when r> 0, 


== 00. when y = 0, ‘ (4.4) 


[ F(r, O)r* dr = A | 
where A is finite, then the integrals can be evaluated giving 


“u” = — 4.5 
48(2r rede vi) — 


f(r, t) ‘in gerere ion ee : (4.6) 
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Loitsiansky describes these initial conditions as referring to a “‘point'source of strength 
A” since the analogous problem in the five-dimensional field is simply the spread of heat 
from an initial point source. Millionshtchikov has also obtained this special solution 
and remarks that it describes the turbulence which exists subsequent to an initial ran- 
dom distribution of concentrated line eddies, provided that the effect of triple correla- 
tions is ignored. Since the triple velocity correlations cannot be neglected in the early 
stages of the decay when w’ is not small compared with the characteristic velocity of 
the turbulence-producing device, Millionshtchikov’s interpretation of the initial condi- 
tions (4.4) cannot be regarded as having physical reality. 

v. Karman and Howarth considered the particular set of solutions of (4.1) which 
are functions of r/(vt)'” only, i.e. a self-preserving solution was assumed. The solution 
(4.6) is the only solution of this kind if certain conditions concerning the behaviour of 
f(r) are accepted. v. Karman and Howarth gave a family of self-preserving solutions of 
(4.1) with the quantity a, = vt/d*, as parameter, viz. 


f(z) -_ ee OM ic. aanke PO. (4.7) 


where x = r/(vt)'””, and M,,,,(z) is the same solution of the confluent hypergeometric 
equation as that defined by Whittaker and Watson (Modern analysis, 1927, p. 337) 
and denoted by this symbol. The solution (4.6) corresponds to the particular value 
a = 4. Now when a < }, the expression (4.7) is proportional to x-*°* when x is large. 
If the restriction that the various moments of f(r) should all be finite is accepted on 
intuitive physical grounds, then certain values of a can be rejected. In particular, if 
the fourth moment {> r* f dr (which occurs in the invariant of Sec. 2) is required to be 
finite then all values of a less than + must be rejected. 

On the other hand, when a > }, the expression (4.7) becomes negative for certain 
values of r since it has the expansion 


_ -x?/8 = (2.5 — 10a)(3.5 — 10a) --- (2.5 — 10a +n — 1) x® * 
I) =e { ba » n! (2.5)(3.5) --- (2.5 +n — 1) (xy (4.8) 





valid for all finite values of x (Whittaker and Watson, p. 337). (Note that although 
the fourth moment fo 7r* f dr converges when a > }, the invariant relation of section 
2 cannot be employed to obtain a definite value for a and thence a definite decay law, 
because the solution makes the fourth moment vanish and A is constant whatever the 
energy decay law). Negative values of f have never been measured, so that there are 
some physical grounds for rejecting values of a greater than }. Thus when these restric- 
tions are applied, v. Karman and Howarth’s family of self-preserving solutions reduces 
to the single solution (4.6). Both restrictions are of course implicit in Loitsiansky’s 
analogy with the propagation of heat in a spherically-symmetrical five-dimensional field. 

Let us return now to the general solution (4.2), with its associated decay law (4.3), 
which satisfies the Eq. (4.1) obtained by neglecting the triple correlation term. When 
tis so large that vt > s” for all values of s for which F(s, 0) is large enough to contribute 
to the integrals, the integrals in (4.3) and (4.2) simplify. Thus when t—> ~, 


" 1 f 4 
2 _, —_____ F(s, 0)s* ds 
. 48(2r)'* (vt)? 0 , 0) 


Se rt 


se ER 
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and 


72 
f] 


1 i . ~r*/Srt 4 
— ————,,, F(s, 0) e s ds, 
48(2m)'(vt)”? 0 , 0) 


. ‘i —r?/Syt 
1.€., f—e Pi 


These limiting forms show that if the correlation f(r, 4) and turbulence intensity u’” are 
calculated on the assumption that the triple correlation is without effect, the solutions 
obtained tend to the forms (4.5) and (4.6) as t > © whatever the choice of conditions 
at ¢ = 0. The common limiting form is a self-preserving solution, so that there is here 
further support for the contention of the previous section that a self-preserving solution 
does exist when ¢ is large, and that it is given by Eqs. (4.5) and (4.6). 

In a later work® Millionshtchikov made an attempt to determine the effect of triple 
correlations at decay times which are not large, using a reiterative method. However 
his attempt is based upon the existence of the solution (4.6) as a first approximation 
and leads only to solutions which have the same self-preserving character. Since a self- 
preserving solution has been shown to be possible only when ¢ is large, it is questionable 
whether Millionshtchikov’s approximate solutions at decay times which are not large 
have any significance.* 

5. Partially self-preserving solutions at Reynolds number which are not large. It 
has been seen that a correlation function which is completely self-preserving can only 
occur when the decay time is large. On the other hand, it is known that some measure 
of self-preservation exists when ?¢ is not large. The function f is always parabolic near 
r = 0. Recent experiments’ have indicated that the expansion of f in powers of r/\ as 
far as the term of fourth degree, and of k as far as the term of third degree, are inde- 
pendent of ¢ at decay times which are not large. This suggests that we should explore 
the consequences of assuming partially self-preserving solutions for the correlation 
functions. We therefore write 

, 


J=IV), hW=kKW, w=s, 


for a range of values of r, 0 < r < l, where 1 is an unknown length. From the above 
evidence / must be at least as great as the maximum value of r for which a fourth degree 
polynominal gives a good representation of f(r). The fundamental Eq. (1.3) again reduces 
to the form (3.2) for the restricted range of r, viz. 


: 4. 5 = 1? dR 1 
(sr + yr 4 By + 59) 53 SP us) + Zee + £8) = 0. (3.2) 


Since the hypothesis leaves arbitrary the behaviour of the correlation functions at 
large values of r, we cannot make use of the invariant relation of Sec. 2. The energy 


°M. Millionshtchikov, On the theory of homogeneous isotropic turbulence, C. R. Acad. Sci. U. R. S. S. 
32, 615-621 (1941). 

*A further criticism of Millionshtchikov’s work is that he has omitted to take into consideration 
the correlation between two velocity components and the pressure when determining the relation between 
triple and quadruple velocity correlations. His idea of determining quadruple velocity correlations from 
the approximation that the relation between double and quadruple correlations is as it is for a Gaussian 
distribution of the sums of velocity components at two points seems to be useful and to warrant further 


exploitation. 
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decay laws are therefore not prescribed as they were in the case of completely self- 
preserving solutions and there will be a degree of indeterminacy in the deductions. 

It was mentioned in Sec. 3 that there are in general two methods of satisfying Eq. 
(3.2) for non-zero ranges of y and ¢. According to the first method, in which the energy 
decay follows a power law, R, is constant and for non-zero values of this constant we 
have n = 1, Le. 

_ 10 


at (5.1) 


u’*? = Bi, d= 10rt, Ry 
where B is a constant. The functions f(y) and k(y) are in this case connected by the 
equation 

gr + {b+ 8) + 5 + hele + 4x) =0 (5.2) 
provided 0 < r < I. The relations (5.1) and (5.2) become identical with Dryden’s de- 
ductions’”’ from the postulate of self-preservation when his scale factor L is replaced 
by the length used in the present analysis, viz. 4; Dryden took Eq. (5.2) to be valid for 
all values of r but we have already seen that the assumption of complete self-preservation 
leads to a quite different set of results. 

It does not seem possible to check Eq. (5.2) relating f and k, since measurements 
of k are difficult to make and no results have yet been published. However, measure- 
ments of wu’ at different stages of decay have frequently been made, and the validity 
of (5.1) can be assessed. The evidence for the law of energy decay has been discussed 
by Dryden.” The data from different sources are not wholly consistent, but inasmuch as 
any one law does describe the experimental relations, it is u’~* ~ ¢", where n lies between 
one and two. In more recent experiments at the Cavendish Laboratory, Cambridge,° 
the decay of \ has been measured simultaneously with that of w’ and the energy equation 
has in this case provided a check on the consistency of the two sets of measurements. 
It was in fact found that the relations (5.1) were obeyed to a quite high degree of ac- 
curacy for the decay range 40 M/U < t < 120 M/U, where M and U are the periodic 
length and velocity associated with the turbulence-producing grid and were varied over 
the range 5.5 X 10° < UM/v < 2.2 X 10*. Under these same conditions it was found, as 
mentioned above, that f(r) and k(r) were self-preserving during decay at least as far 
as the terms in r* and r* respectively. The above experiments are therefore quite con- 
sistent with the postulate of limited self-preservation and with the deductions obtained 
by using the first method of satisfying Eq. (3.2). A point left open is the value of 1, 
which presumably can be determined by the region of validity of the relation (5.2). 

According to the second method of satisfying Eq. (3.2) (see Sedov’), the coefficients 
of bracketed terms therein are related by 

2 
a oa = aR, + b (5.3) 
where a and 6 are constants. For solutions of this equation which do not make R, con- 
stant, f and k must then satisfy the equations 


“HH. Dryden, [sotropic turbulence in theory and experiment, v. Karman Anniversary Volume on Ap- 
plied Mechanics, 1941, pp. 85-102. 
1H. Dryden, A review of the statistical theory of turbulence, Q. Appl. Math. 1, 7-42 (1943). 
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4 
k’ + Tk = ~ayf’, (5.5) 


provided 0 < r < lL. Apart from the constants a and b, this method of solution provides 
unique determinations of the correlations f and k within the restricted range of r.* 
Consider the decay law described by (5.3). This equation is to be solved with the 
aid of the energy equation in the form 
2n° dR, dy? /y 


»R, Y = “at — 10. (5.6) 


The solution cannot be obtained explicitly, and is given by 
2 2 2 
A Ry vR) 


- = K(a + b/R,)*°*”” = u” (5.7) 





where K is a constant of integration, and 


dR, oi ( »ye" 
a ~ Ket (5.8) 


It will be seen in a moment that a is necessarily negative and it is then evident from 
(5.7) that b must be positive if u’ is to become large for some values of R,. The nature 
of the decay relations specified by (5.3) is as indicated in Fig. 1. The variation of wu’ 


~bh... a ee sons 








Fie. 1. 


has the general features of measured energy decay curves. R, increases with decay time 
and asymptotes to a constant value —b/a. It follows that u’~* and \* asymptote to the 
values (—a/b) 10t/v and 10 vt respectively as ¢ increases. 





*This might be regarded as a hint that the solution is physically impossible. However I have not 
been able to find any definite anomalies. 
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It might be argued that this method of satisfying Eq. (3.2) has thus far led to pre- 
dictions which are not inconsistent with experiment, since we cannot be sure that the 
asymptotic variation (2, = constant) does not occur over the range of decay times 
used in the experiments mentioned earlier. But the Eqs. (5.4) and (5.5) determining 
the functions f(y) and k(~) make possible a further comparison with experiment. There 
is available sufficient experimental evidence to determine the value of the constant a@ 
in (5.5). It has been found’ that when r is sufficiently small, the function k(r) has the 


ii «as s( ") (5.9) 


6 A 


form 


where S is an absolute constant (of value about 0.39) for the ranges of decay times 
and mesh Reynolds numbers used in the experiments. There are also theoretical reasons, ’” 
derived from Kolmogoroff’s theory of locally isotropic turbulence,’* for believing that 
S is an absolute constant whenever the Reynolds number is sufficiently high. Thus 
comparison of the coefficients of powers of r* in (5.5) shows 


a 


a=-—-S. (5.10) 
6 
Then if R{ be written for the asymptotic value of the Reynolds number Rh, , = u’d/», 
the value of b is 
b = —ah{ = - SRi. (5.11) 
) 


Changing the variable of (5.4) to x, = ¥/(a)'”’, where 


l 
“10+ 2b 10 + 78Ri/3’ 
the equation for f becomes 
{ 
fv 4 r( 4 x) + Saf = 0 (5.12) 


v. Kaérmadn and Howarth’ have pointed out, in a slightly different context, that the 
solution of this equation is related to the hypergeometric function and that the solution 


which satisfies f = 1 and f’ = 0 when r = 0, can be written 
*/ I'(5/2) [ 10a—1 3/2—10a_ ~—(x*r)/8 
f(x) = el = 2) oe 5.13 
IM '(10e)1'(5/2 — 10a) J, * ( 7) , — (5.13) 


when, as is here the case, a < }. Knowing the value to which u’d/v tends‘according to 
this method of satisfying the requirement of limited self-preservation of the correlation 
functions, it is thus possible to determine f. Figure 2 shows not f, but the related double 


correlation function g, which, in the notation of Sec. 1, is given by 
ul’ g(r) = (uyus),-; - 


€i=0 


143. K. Batchelor, Kolmogoroff’s theory of locally isotropic turbulence, Proc. Camb. Phil. Soc. (to 
be published 

8A. N. Kolmogoroff, The local structure of turbulence in an incompressible viscous fluid for very large 
Reynolds numbers, C. R. Acad. Sci. U. R. S. 8. 30, 301-5 (1941) and 32, 16-18 (1941). 
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g(r) is easier to measure in practice, and is related to f(r) by the continuity equation 


r af(r) ; 


ar) = fr) +55, (5.14) 


g has been calculated from (5.13) and (5.14), in part by numerical integration and in 
part by the use of an asymptotic expansion of the integral, for a number of values of 
R{ comparable with those found in the experiments at Cambridge. The largest value, 


viz. 44, is the value of u’\/v measured at intermediate decay times with a mesh of circular 
rods at 1’’ spacing and a wind speed of 42 ft/sec. 
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Fig. 2. Solutions of d?f/dy* + (4/y + 4(5 + b)y) df/dy + 5f = Oandg =f + jy df/dy, 
where y = r/d, b = 7SRX/6, S = 0.39. 


Mr. A. A. Townsend and the author hope that it will soon be possible to present 
measurements of g(r) during decay under different conditions in order to determine the 
validity of the solution (5.13). Assuming that the measurements show u’\/v to be con- 
stant at decay times which are not too large and that self-preservation of the function 
f(r) during decay occurs over a limited range 0 < r < I, then agreement between ex- 
periment and the curves of Fig. 2 over this same range* of values of r would show that 


*A hint concerning the possible range of validity of the curves in Fig. 2 is provided by the require- 
ment, deduced from (5.14), that the complete function g(r) obeys the relation 


is rg(r) dr = 0. 
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the second method of satisfying Eq. (3.2) is here valid. Necessary consequences are 
that k is given by (5.5), i.e. 


ou fe 
k(y) = . Sy / y 5 dy O<r<J (5.15) 


where f is given by the expression (5.13), and that at decay times which are sufficiently 
small the decay of wu’ is such that u’d/v increases from a small value to its asymptotic 
value Ri . On the other hand, lack of confirmation of the curves of Fig. 2 would show 
that the first method of satisfying Eq. (3.2) is valid, i.e. that u’d/v is constant over the 
whole of the decay range for which limited self-preservation holds and that k is de- 
termined by (5.2) as 


“y) = —2 gy-* [ ys at 
Rk(y) = -2 af + vf) — 3y | vay (5.16) 


The function f(y) is not prescribed by this method of satisfying (3.2), except in the 
particular case mentioned below. 

The case in which R, , or Ri in the second method, is very small (but constant) is 
particularly interesting, for the two methods of satisfying (3.2) then lead to identical 


approximate equations for f(y), viz. 
4,5 
i +9(5 + 5¢) + 5f = 0. (5.17) 


The solution is shown in Fig. 2 and it can be predicted that the correlation function 
f(r) will approximate to this shape at low Reynolds numbers of turbulence over the 
range of values of r for which this same correlation preserves its form during decay. 

6. Solutions at high Reynolds number. (a) No assumption of self-preservation. Using 
the energy Eq. (1.4), the basic Eq. (1.3) can be written 


*) of (10f af’) 
, my Be actad = — canes — a? ~ : 
u (i + ig ry W>2 +f’ + (6.1) 
We compare the orders of magnitudes of the two terms 

10f and f+ il 

A ; r 


When r is sufficiently small for f to be represented by the parabola 1 — r°/2)’, these 
terms are of equal order of magnitude, however small \ may be. Now when the Reynolds 
number of the turbulence is increased observation shows that the correlation curve 
does not change appreciably outside the parabolic region, the change being confined to 
a diminution in the value of \. The first of the two terms will therefore become relatively 
large at high Reynolds numbers for values of r lying beyond the parabolic region, and 
the approximate form of (6.1) is 
ak | of  10rf 

, , . ? a] 
u\k’ + - =< — —~. 6.2 
( r ot Ss 8.3) 
It is possible to integrate Eq. (6.2) if r is further restricted to be sufficiently small 


for the approximation f = 1 to be valid. In this case, 
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4k —-:10v 10 6.3) 


ee - a 
and the solution which makes k(r) vanish at the origin is 
k(r) = -> >. (6.4) 


The range in which (6.4) can be expected to apply is rather restricted. In the first place, 
the Reynolds number must be large enough for (6.2) to apply. r must be large enough 
to lie outside the parabolic region of f(r), but in view of the requirement of high Rey- 
nolds number this is not an important limitation. More important is the requirement 
that r should be small enough for the value of f to be close to unity. The meaning of 
(6.4) is evidently that at large Reynolds numbers, the region near the origin in which 
k(r) is cubic becomes small and the curve tends to a straight line with a slope which is 
determined by the turbulence Reynolds number R,. This discussion is substantially 
that already given by Kolmogoroff.’* 

(b) Self-preserving solutions. Under the assumed conditions of high Reynolds num- 
ber, the dissipation length parameter \ is small compared with other lengths associated 
with the turbulence. In the limit, \ is zero and the tangent to the correlation function 
f(r) at r = O is not horizontal, but makes an angle with the abscissa which is probably 
90°. In discussing self-preserving solutions at high Reynolds numbers it is no longer 
possible to use A as the reference length and some other length LZ associated with f(r) 
will be employed. We specifically permit \/Z to vary during decay, since we should 
otherwise merely repeat the analysis of Sec. 3. 

The present hypothesis is thus 


f) = $0), KO =k, x=F (6.5) 


for a range of values of r to be specified later. The reference length Z must of course 
be defined by values of the function f(r) within the range of r for which (6.5) is to hold. 


The basic equation (1.3) becomes 


LdL df 1 (2 4k) (4 4 af) 
— aie cadeaiae ——_—_ =— — ania } mies } —, } =n aaiie 
aa Yak Qv dt * dx 2 R, dx x dx” x dx (6.6) 
where R, = u’L/v. 


L/ is large under the assumed Conditions, and the term (d’f/dx? + 4/x df/dx) 
should be neglected by comparison with the term 5fL*/\* in order to be consistent 
with the assumption that the non-self-preserving parabolic region near the origin is 
infinitesimal in extent. Hence (6.6) becomes 

L? L dL ( at) 
Pal a ea , R,\ k’ —]} = 0. 6.7 
a) + Se Oh) + Rak’ +o 6.7) 

Equation (6.7) will be satisfied if the three coefficients, which are functions only of 

t, are proportional, i.e. if 


Ri, = AL’/x, (6.8) 


L dl 
vy at 


= BR,. (6.9) 
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These are the equations derived by v. K4rman’ in a discussion of self-preserving solu- 
tions at high Reynolds number. The equation connecting f and k becomes 


10f + ABxf’ + Ai + : k) ~% (6.10) 


As before, it is also possible to satisfy Eq. (6.7) in the manner suggested by Sedov,’ but 
in this case there are reasons for rejecting this alternative. The appropriate relation 
between the coefficients is 

L? * Ld 
10> = a——= + dR., 11 
Pp "*at oy 
where a and b are constants. For solutions of this equation which do not make L dL /v dt 
proportional to R; , the functions f and k must then satisfy 


af + xf’ = 0, (6.12) 


bf +k + : k= 0. (6.13) 


However Eq. (6.12) has no solutions such that f = 1 when y = 0 and the method of 
satisfying (6.7) must be rejected. 

v. Kdrmén’s Eqs. (6.8) and (6.9), and Eq. (6.10), therefore appear as necessary 
consequences of self-preservation at high Reynolds numbers. The meaning of (6.8) 
becomes clear when we substitute for \” in the energy equation to obtain 


This equation shows that the quantities wu’ and L may be considered as characteristic 
of the whole turbulence in a calculation of the work done against Reynolds stresses, as 
indeed is to be expected when the correlation functions preserve their shape during 
decay. Perhaps less to be expected is that this expression for the energy decay remains 
valid (as will be seen later) when the correlation functions are only partially self-pre- 
serving. 

When + is small and f may be replaced by unity, Eq. (6.10) becomes 


dk , 4 10 °° Ll’ 
dx +" nr as 10 aR 





which is identical with Eq. (6.3) derived without any assumption of self-preservation. 
Clearly the general solution (6.4), viz. 
2r 
kr) = ——-, 

”) Ryd 
is self-preserving with L as the unit of length when the decay law (6.8) holds. 

The decay equations (6.8) and (6.9) contain three variables, Z, A, u’ and may be 

solved with the help of the energy equation (1.4). The differential equation for L is 


CL 5 (2) 
L d?® AB gh (6.14) 
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and has the solution* 


t (AB)/(5+AB) 
L = L,(“) ° (6.15) 
Then from (6.8) and (6.9) 
£ * (5 4. AB)t, (4)" 5+AB) si 1 — 
u’ ng A | . — u; ti ’ (6.16) 
\” = o(5 + AB)t, (6.17) 


where LZ, and uj, are the values of L and w’ att = t . 

It is at this stage that we must be more specific about the range of values of r over 
which self-preservation is to be postulated. If the correlation functions are completely 
self-preserving, Loitsiansky’s invariant relation (2.1) shows that 

u”L> = constant, (6.18) 


and we must have 


AB = 2, 


as has been pointed out by Kolmogoroff.'* The decay laws then become 


Ea (8)" 4. (t)” fs 
Lo a VA ‘ u’ as ty 4 A . Tt. (6.19) 


The corresponding relation between f and k is, from (6.10), 
: F q 
10f + 2xf’ + Ake + . k) = 0 (6.20) 


which can be integrated to give 


pa - 5 y= - (2% ) rs (6.21) 


Titus L 


At the present time, we have not got sufficient evidence at high Reynolds numbers of 
turbulence to determine the validity of the decay laws (6.19) and the correlation relation 
(6.21). 

If, on the other hand, the hypothesis of partial self-preservation only is made, the 
value of AB remains arbitrary and deductions about the decay laws finish at the rela- 
tions (6.15)—(6.17). We can consider one or two consequences of particular values of 
AB. The assumption of large Reynolds numbers entered the analysis via the postulate 
that L/d is large. Since (6.15) and (6.17) show that 

4 AB—5)/(2(5+AB)) 
+ = is (t) (6.22) 


d [5 + AB)vt.]'”? \to 





the subsequent behaviour of L/\ depends critically on the sign of AB — 5. 


*y. KA4rmAn! has written the exponent of the solution in error as 5/(5 + AB). 
4A. N. Kolmogoroff, On degeneration of isotropic turbulence in an incompressible viscous liquid, 
C. R. Acad. Sci. U. R. S. 8. 31, 538-540 (1941). 
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If AB > 5, L/X (and also w’L/v) increases and the assumed state of affairs applies 
with ever-increasing accuracy. A state of turbulence in which the Reynolds number 
u'L/v increases indefinitely—although the energy decreases—does not seem possible, 
and an argument of the kind used in Sec. 3 does in fact confirm this impression. Equation 
(2.3) can be written 


9 fan f™r ocr Bul? > or (6u” | 
: (ns | —fd- = =~, | kd ; i ——, (6.23) 


at Vv 0 v 


Suppose that when t — ©, we may write 


i L ds L mrs fg k dF ~ 0 


and the decay laws (6.17)-(6.19) are valid. Replacing the terms of equation (6.23) by 


their orders of magnitude in ¢, 


o(t'?4" 10)/(5+AB)+a * rea of** eres 4. _ iamasiaaiiaamasas * 


In the present analysis the viscous term (i.e. the last on the right side) is assumed to 
be without effect, so that the remaining terms must be of equal order, i.e. a = 6. More- 
over, with the assumption made previously that {> k dr < 0, the term on the left of 
(6.23) must be negative and consequently 


10 — 2AB 


5 + AB’ 


The first term on the right side of (6.23) thus varies as a power of ¢ which is less than 
(—1), and the viscous term will only be of smaller order if AB < 5. Since L/d increases 
indefinitely with t when AB > 5, we have here an inconsistency which prohibits such 
values of AB except at small values of t. 

The case AB = 5 leads to the decay laws 


ul’? ~t, L’ ~ »* = 10 vt. 


Since L/X is constant, there is no tendency for the approximation on which the solution 
is based to become invalid as ¢ increases. Nevertheless such decay laws cannot persist 
indefinitely. When ¢ is large, the first term on the right side of (6.23) must vary as some 
power of ¢ not less than (—1) if the viscous term is not to become dominant; however 
such a variation with ¢ requires [¢ (r/L)fd(r/L) to approach — ~, which is impossible. 
Finally, when AB < 5, L/X diminishes as ¢ increases so that there comes a time when 
the postulate of high Reynolds number ceases to be valid; the case AB = 2 deduced 
from complete self-preservation is in this category. The system is here naturally unstable, 
whereas if the case AB > 5 occurs at all some secondary factor disturbs the system 


before ¢ becomes large. Thus none of the regimes deduced for large Reynolds number 


are possible for indefinitely large times of decay. 
7. Conclusion. The various hypotheses discussed in the last three sections, and the 
deductions made from them, are summarized in the following table. 
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comple te self- 
preservation of 
f(r) and k(r) 


during decay 























self-preservation 
of f(r) and k(r) 
for0O < r<l 
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Reynolds number not large Large Reynolds number 
af 
2 wien < 
rs oe ; __ for r so small that f(r) 
w?~ 7, -d? = vt, 
fwW=e" ns 
r 
— where y = x 
went" Dot”, = Tt 
Only possible when t ©; see 
entry above. (7) =— 2 ( Ze | dl (5) 
L 7 \tuz/ L°\L 
‘Either a a yl? ww p10 54 AB) 2, pPAB/(5+ AB) 
(1) u’* ~t, dr? = 10pt \” = (5 + AByt 
f" + Ae +: SY) + 5f + 
+ Ak +4) =0 
+1 pe + +2 ti) = 0 x 
O<r<) where x = 7, 0O<r<l. 
or 
b 10/b 
(2) = Kila + | . 
- 1 “ a 
t= K | n — R, dR, 
” , $ ¥) 
prts(ie 
5 
cxize-* 
10 ve 0.91 Rx ° 
0.46 _ 
Ey) = 948 [” ype ay, 
y 0 
where Rf = — 2 G<r <€ ¢ 
In either case, 
4 
preset 5H) 4 sy 
= 0 O<r<J 
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In conclusion, it should be noted that even if the deductions listed in the table are 
found to be in agreement with the appropriate measurements, the analysis of this paper 
does not by any means constitute a solution to the problem of the decay of isotropic 
turbulence. ‘All that can be inferred from agreement with experiment is that the decay 
process takes place in a certain manner, for example, with self-preservation of the 
correlation functions for a certain range of r and all that this implies. (The deductions 
about the turbulence as t — ©, or when r is small and the Reynolds number is large, 
are exceptions to’ this statement since no hypotheses were required in these cases). 
In other words, the hypotheses of the above table are mathematical in origin and, so 
far, have no physical raison d’étre. The task now is to find physical reasons why the 
decay process for ¢ not large takes place as it does; indeed this might be said to have 
always been the chief problem of research on isotropic turbulence. Kolmogoroff’s simi- 


. 13,12 . : . . 
larity hypotheses are physical notions which promise to “ 


explain”? very success- 
fully the structure of the turbulence (at high Reynolds number) at any instant. When 
we obtain equally convincing physical ideas about the way in which the governing 
parameters of this structure change from instant to instant, then our progress in the 
problem of isotropic turbulence will be considerable. 
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THE GENERAL PROBLEM OF ANTENNA RADIATION AND THE 
FUNDAMENTAL INTEGRAL EQUATION, WITH APPLICATION TO 
AN ANTENNA OF REVOLUTION—PART I* 


BY 
G. E. ALBERT 
Naval Ordnance Plant, Indianapolis 
AND 
J. L. SYNGE 
Carnegie Institute of Technology 


Introduction. The main result of Part I is the integral equation (3.17)—a general 
and exact equation satisfied by the current flowing in a perfectly conducting radiating 
antenna, the surface of which is a surface of revolution, but otherwise unrestricted 
except for implicit conditions of smoothness. Part II deals with the solution of that 
integral equation. 

No use is made of vector or scalar potentials; an approach through the electro- 
magnetic vectors appears simpler and more direct, since a radiating field is determined 
by the values of the tangential component of the electric vector on the inner boundary 
of the infinite region.’ 

Section 1 gives the basic equations and boundary conditions, together with some 
formulae connected with a dipole. Section 2 deals with the case of radiation in a finite 
cavity; this simpler case points the way to the theory of antenna radiation. The basic 
integral equation (2.6) for cavity eigen-modes is obtained, as well as the integral equa- 
tion (2.4) which has the same form as the basic equation for the exterior (antenna) 
case. Section 3 develops the general integral equation (3.1) or (3.2) for antenna radiation, 
which reduces to the comparatively simple explicit form (3.17) for an antenna of revolu- 
tion.” In Sec. 4 some actual radiating systems are discussed from the standpoint that 
the tangential component of the electric vector determines the field. 

Section 5 is devoted to a discussion of the gap. The gap is in a sense the most im- 
portant feature of an antenna because it is the source of radiant energy. (No energy 
passes out through the surface of a perfect conductor.) Strangely enough, the gap has 
received scant attention in radiation theory. The recent work of Infeld* shows that 
the question of the gap is a most delicate one; a very narrow gap leads to zero impedance 
unless there is a compensating thinness of the antenna at the gap. Zero impedance is 
avoided in the present paper by keeping the gap long compared with the thickness of 
the antenna at the gap as in (5.7). 

In Sec. 6 the unknown is changed from the current /(z) to the relative current $(z) , 


*Received March 31, 1947. This paper is based on a report written by one of the authors (J.L.S.), 
when at the University of Toronto. This report was issued in 1942 by the National Research Council of 
Canada. The method now used in Part I differs considerably from that used in the report. 

1The importance of the tangential components of the electromagnetic vectors appears to have been 
first pointed out by H. M. MacDonald, Electric waves, Cambridge, 1902, p. 16. 

2For other, and we believe less exact integral equations in the antenna problem, see E. Hallén, Nova 
Acta Reg. Soc. Sci. Upsaliensis, 11, No. 4 (1939); P. Nicolas, L’Onde électrique, 18, 193-211 (1939); 
L. Brillouin, Quarterly Appl. Math. 1, 201-214 (1948). 

3L.. Infeld, Quarterly Appl. Math. 5, 113-132 (1947). 
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viz. the ratio of the current at z to the current at the center of the gap. Formulae (6.5), 
(6.14) give the impedance Z in terms of ¢(z). The function ¢(z) satisfies the integral 
equation (6.3), not essentially different from (3.17). 

The thin antenna is discussed in Sec. 7. We are able to show that the current is 
approximately sinusoidal outside the gap (7.10). The sinusoidal distribution of current 
is commonly accepted, but proofs up to the present have not been satisfactory. The 
present method is valid for any thin antenna of revolution; it is not assumed to be 
cylindrical, except at the gap. 

The calculations in the case of a finite gap are too complicated to carry out at present, 
and for simplicity the gap is assumed to be short in Sec. 8. Whatever the practical 
validity of the final results may be, it must be clearly understood that for theoretical 


validity we must have 


wave length > length of gap > radius of antenna. 


Equation (8.10) gives in compact form the principal part of the impedance of any thin 

antenna with such a gap anywhere in it. If the gap is taken at the center and the shape 

term X, neglected—it is negligible for a cylindrical antenna with plane ends—we get 
the formula obtained by the less rational method of Labus.* 

In Sec. 9 the shape term is discussed, with reference to the following thin antennae: 
(i) a cylinder with spheroidal ends, (ii) a spheroid, (iii) a cone. It is remarkable that 
the half-wave spheroid has zero reactance. 

Section 10 contains a method of successive approximations for dealing with any 
antenna, thin or thick. The first approximations for the current and impedance agree 
with the principal parts obtained in Sec. 7 for the thin antenna. Higher approximations 
have not been calculated. 

The principal features of the paper may be summarized as follows: 

1) Radiation problems are reduced to the solution of integral equations involving the 
tangential components of the electric and magnetic vectors. 

2) The gap is treated more adequately than has been done before. 

3) The approximately sinusoidal current distribution for a thin antenna is established. 
The usually troublesome question of end effects disappears automatically in the 
present method. 

4) The usual first approximation for the impedance of a thin antenna is obtained with 
the following generalizations: . 

(a) The antenna is not necessarily cylindrical with plane ends. 

(b) The gap is not necessarily at the center. 

5) A method of successive approximations is given for dealing with thick antennae. 
Advice and criticism from the following is gratefully acknowledged: Professors L. 

Infeld, V. G. Smith, and A. F. Stevenson, and Dr. W. Z. Chien, at the University of 

Toronto; and Professors A. E. Heins, and A. Weinstein, at the Carnegie Institute of 

Technology. Computations for the graphs were made by Mr. George Succop. 

1. Notation and preliminary equations. We confine our attention to simple harmonic 
electromagnetic fields in vacuo; the electric and magnetic vectors are the real parts of 


E exp (—iket), H exp (—vket), (1.1) 


4J. Labus, Hochfrequenztechnik und Elektroakustik, 41, 17-23 (1933); G. H. Brown and R. King, 
Proc. I. R. E., 22, 457-480 (1934). 
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where E and H are vector functions of position, in general complex, c is the velocity of 
light, and & is a real constant related to the wave length \ by \k = 2z. Except for the 
statement of some final results, we shall use Heaviside units, so that Maxwell’s equations 
read 

—ikE = V Xi, ikH = V XE, (1.2) 
V-E = 0, V-H = 0. (1.3) 

Equations (1.3) are necessarily satisfied if (1.2) are satisfied. 
As regards boundary conditions, we have on the surface of a perfect conductor, 


with unit normal n, 


nXE= 40, (1.4) 


i.e. E, = 0, where E, is the tangential vector component of E. It follows from the second 
of (1.2) that on the surface of a perfect conductor 


H-n = 0. (1.5) 


In dealing with an exterior problem, conditions at infinity must be stated. Instead 
of the usual Sommerfeld conditions, we shall take the following as conditions for out- 
ward radiation at infinity. A constant B exists such that at all sufficiently great distances 
p from some fixed point O we have 


|E| < B/p, |H| < B/p, (1.6) 
|E+nxXH| < B/p’, |\H—nxXE| < B/p’, (1.7) 


where n is the unit vector drawn from O towards the point at which the field is considered. 
It is clear that these conditions are invariant under bounded displacement of the 
point O; this causes a bounded change in p and a change in n of order 1/p. 
The following extended Poynting’s theorem, due to H. A. Lorentz,’ is basic in our 
work: Let (E, H) and (E’, H’) be two Maxwellian fields in a finite region R bounded by a 
surface S, and let neither field have any singularities in R or on S. Then 


| n-(E X H’ — E’ x H)dS = 0, (1.8) 
Ss 


where n is the unit normal vector to S. This is proved very easily by Green’s theorem 
and Maxwell’s equations (1.2). 

Consider now the case where the region R extends to infinity, and is bounded intern- 
ally by one or more surfaces S. Then if = is any large sphere of radius p enclosing S, 
we have 


[o@xH-E xmas=/[n&xH’- Xd, (19) 


“§ « 


where n denotes the outward unit normal vectors to S and to ¥. Suppose that both 
fields satisfy the conditions (1.6), (1.7) for outward radiation at ififinity. Then, denoting 
by 7 terms of order 1/p*, we have on 2 





‘Cf. Frank and Mises, Differential- und Integralgleichungen der Mechanik und Physik, Vol. 2, 


Braunschweig, 1927, p. 576. 
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n-(E X H’ — E’ X H) = —n-[(a X H) X H’ — a X OH’) X H+ T; 


—n:(H(n-H’) — n(H-H’) 
(1.10) 


— H’a-H) + n(’-H)) + T; 


= Zs . 
Hence, as p tends to infinity, the integral over = tends to zero, and we have this result: 
Let (E, H) and (E’, H’) be two Maxwellian fields in an infinite region R bounded internally 
by one or more surfaces S; let neither field have any singularities in R or on S, and let each 
field satisfy the conditions (1.6), (1.7) for outward radiation at infinity; then 


[ .-@x H’ - Bx Has =0. (1.11) 
8 


Suppose now that there is a dipole of vector strength A at the point P’; the field 
due to it at a point P is 

E’ = V(A-Vyv) + RYVA, H’ = ikA X Vv, y = PP’ / PP". (1.12) 

This field plays an important part in our work. Its essential features are: (i) it satisfies 

Maxwell’s equations, (ii) it has a pole of a certain type at P’, (iii) it satisfies the condi- 

tions (1.6), (1.7) at infinity. Let (E, H) be any Maxwellian field, regular in the neighbor- 

hood of P’, and let o be a spherical surface with center at P’ and small radius p. Then 


lim | n-(E X H’ — E’ X H) dS = 4mikE-A, (1.13) 
p-0 ¢ 
where E on the right hand side is evaluated at P’, and n is the unit outward normal 
to «. This limit was incorrectly calculated by Sommerfeld;* the error has been pointed 
out by A. F. Stevenson in a paper in course of publication,* and the correct limit (1.13) 
given. 
2. The integral equation for cavity radiation. Interior (cavity) and exterior (antenna) 
problems in electromagnetic radiation present some common features and some striking 


ae p’ e 
R 


Ss 
Ss 
Fic. la. Interior Problem: Dipole Fig. lb. Exterior Problem: Dipole 
at P’ yields integral equation. Dipole at P’ yields integral equation. Dipole 
at P” yields field at P”’. at P” yields field at P”’. 


6P. Frank and R. v. Mises, op. cit., 578-580. 
*Footnote added in proof: A. F. Stevenson, Q. Appl. Math. 5, 369-384 (1948), in particular p. 377. 
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differences. Although the problem with which we are primarily concerned is the problem 
of antenna radiation, we shall first discuss cavity radiation. 

Let R be a finite region (cavity) bounded by a surface S (Fig. la). Let (E, H) be 
a Maxwellian field inside R, satisfying (1.2); so far we do not mention boundary condi- 
tions. Let (E’, H’) be the field due to a dipole situated at P’, outside R. Then, as in (1.8), 


[ n-(& x H’ — Fx B) dS = 0. (2.1) 


“S 


Let us write 


E,=nXE, H,=2xX4H, (2.2) 


so that E, , H, are vectors tangent to S; they are respectively the tangential vector 
components E, , H, of E, H, turned through a right angle about n in the positive sense. 


Then (2.1) may be written 


| &,-H’ + H,-E) as = 0, (2.3) 


“Ss 
or 


0, (2.4) 


| &,-H) + H,-E) ds 


where E! , H/ are the tangential vector components of the dipole field. 

Noting that the vector strength and the position of the dipole are arbitrary (subject 
only to the condition that it must lie outside R), we see that (2.4) is an integral equation 
satisfied by the tangential components on S of any Maxwellian field regular inside S. 

Once (E, , H,) are known on S, we can find the field (E, H) at any point P” in the 
cavity as follows. We set up a dipole of strength A at P’”’; let (E’’, H’’) be the field of 
this dipole. Then we apply (2.1) to the region bounded outside by S and inside by a 
small sphere centered at P’’. Then, if n is the unit normal to S, drawn into the cavity, 


we have by (1.13) 
4mik(E)p.-A = — | (B,-Hy’ + H,-Ey’) ds. (2.5) 
“S 


The integral equation (2.4) may be used to study forced oscillations of a cavity, 
or free oscillations. Suppose, for example, that the cavity is excited by given E, on its 
surface; then (2.4) becomes an integral equation to determine H, , and, if it can be 
solved, the field in the cavity is given by (2.5). In the case of free oscillations in a cavity 
with perfectly conducting walls, we put E, = 0; then (2.4) reads . 


| H,-E/ as = 0. (2.6) 


It is well known that free oscillations in a cavity are possible only for definite values 
(eigenvalues) of k; (2.6) is an integral equation of an unusual type for the determination 
of such eigenvalues. But we shall not discuss the question of eigenvalues here.’ 


7Cf. H. Weyl, Journal fiir Mathematik, 143, 177-202 (1914); Rendiconti del Circolo Matematico id 
Palermo, 39, 1-50 (1915): 
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The reader should remember that, in the above work, 
(E’, H’) = field due to dipole outside the cavity; 
(E”, H’”’) = field due to dipole placed in the cavity. 


3. The integral equation for antenna radiation. Let R be an infinite region bounded 
internally by a surface S. We are interested in simply harmonic Maxwellian fields, 
satisfying the conditions (1.6), (1.7) of outward radiation at infinity and certain condi- 
tions over the surface S (Fig. 1b). In the physical problem, S is a surface enclosing the 
antenna or coincident with it. 

Let (E, H) be the field under discussion. Let (E’, H’) be the field of a dipole of vector 
strength A situated at any point P’ inside S, i.e. not in the region R. We may apply 
the result (1.8) to the surface composed in part of S and in part of a large sphere. As 
we have already seen, the contribution from the large sphere tends to zero as its radius 
tends to infinity. Hence, in the notation of (2.4) we have 


[ (E,-H! + H,-E/) dS = 0, (3.1) 


where ¢ signifies ‘‘tangentia] vector component” and E, , H, denote that E, , H, are 
turned through a right angle in the positive sense about’ the unit normal n to S. We 
shall in future regard n as drawn out from S, i.e. into R. It is clear that E,-H! = —E,-H? ‘ 
and so (3.1) may be written 


[ H,-E/ as = [ E,-H’ dS. (3.2) 


d 


Equation (3.1) or (3.2) is the fundamental integral equation of antenna theory. If 
this equation is solved, and E, , H, found, the field at any point in & is obtained by 
taking a dipole of vector strength A at any point P” inside R, i.e. outside S. Then, as 


in (2.5), we have 
4rik(E)p -A = — / (E,-H!’ + H,-E/’) dS, (3.3) 
8 
when (E’’, H”’) is the field due to the dipole placed at P’”’; the normal on S points outward. 


Let us now apply our equation (3.2) to the case of an antenna having symmetry of 
revolution about the z-axis (Fig. 2). 














Ss Y Xx 
o7) 
R + 
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© —) > = rs 





Fic. 2. Antenna of revolution with gap. 





1948 THE GENERAL PROBLEM OF ANTENNA RADIATION 123 

The antenna is broken by a gap containing a generator. Let the surface of the antenna 
be extended by a mathematical surface G covering the gap, and let S be a closed surface 
consisting in part of the surface of the antenna and in part of the surface G drawn over 
the gap. We shalJl suppose the excitation of the antenna to be symmetric in the sense 
that, on G, E, lies in the meridian plane. It is natural, then, to consider only fields of 
radiation having this type of symmetry, i.e. E at any point P lies in the meridian plane 
through P and the z-axis, and H in consequence is perpendicular to the meridian plane. 
On S, H, lies in the meridian plane. 

Now take a dipole at any point P’(z.) on the axis of the antenna inside S, the axis 
of the dipole lying along the axis of the antenna, and its strength being unity. We shall 
apply the fundamental equation (3.2). 

Introducing cylindrical coordinates (z, R, w), we have, by (1.12), 


y =e"'/r, r = R? + (z — 2)”. (3.4) 
The surviving components of the field of the dipole are 
» 9, eo » — 3, 2¥ 
Bima thY, Ea=spay He = % ap @.5) 


where the differentiations of y are carried out with R, z, z treated as independent 
variables. However, on the surface of the antenna, RF is a function of z, and so we may 
write Y = y¥(z, 2,). Let us temporarily employ d as a symbol of partial differentiation 
when y and its derivatives are regarded as functions of z and z, . Then, noting that 


ddy_ 9 dpy_ _ 9 dy 
dz,9@R d%0R 2 29R’ 
doay_d9dy_ _ daddy 
dz, 02 02 02 #£O2 dz’ 


we have 


dy odyvdR , dy 
(3.6) 





The tangential component E/ may be written E{ = E’t, where t is a unit tangent vector 
to the meridian curve, as shown in Fig. 2; t makes an angle x with the z-axis, so that 


. ial Sa dR /dz ‘a 1 
m x (1 + (dR /dz)?}"””” OS x fl + (dR/dz)}"™" (3.7) 


Then, by (3.5) and (3.6), 
E’ = E! cosx — Ez sinx 








27-1/2 vy * ay az | 
(1 + (@R/dz)’] cE: + EW YT OR az de 





(3.8) 





(1 + (anjasyy""| BY -£ ¥ | 
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If we write H, = Ht, we have 





— _ 2 ~ 
H,-E; = E’H = [1+ (aR jae)" By - 4 ja. (3.9) 
0 
Also, we may write E, = Ft, H = H’t, where, by (3.5), 
7 — ye = op OW / 
H’ =H. = i aR’ (3.10) 


Now yw is a function of r, a single variable, and dy/dr will denote the derivative. We have 


dR dr dR r dr’ 
and so 


E,-H, = EH’ = ikE Rdy (3.12) 


r dr 
Thus our basic equation (3.2) reads 
[ Ee i oe ht 4+ (dR/dzY}"'?H dS = ik ‘ “ sr E dS. (3.13) 
Let us denote by J(z) the current flowing on the surface of the antenna across the 
section z = constant. Then 
I(z) = 2ncRH. (3.14) 
We shall regard this as the definition of J(z) on G. Also 
dS = 2rR{1 + (dR/dz)*]'” dz. (3.15) 


Hence (3.13) may be written 


/ (iy e fre ds = Selle / RU + (dRfdz)}” t o E(@) dz. (3.16) 

This equation (3.16) is the basic integral equation for the determination of the 
current in an antenna of revolution. We shall therefore restate it in complete form 
with slightly revised notation: Consider a perfectly conducting antenna having the z-axis 
for axis of symmetry, and excited symmetrically. Let the antenna extend from z = |, to 
2 =1,. Let R = R(z) be the equation in cylindrical coordinates of a surface S which consists 
of the surface of the antenna and the gap G. Let E(z) be the tangential component of the 
electric vector on the gap, and E(z) = 0 outside the gap. Let I(z) be the current across the 
section z = const. as given by (3.14). Then I(z) satisfies the integral equation 


ls 
/ K(z, 2.)I(z) dz = iM(z), L<ae<i,,. (3.17) 
bi 
where 
, _ _ ov ps 
K(z, 2) = de Oz + ky, 
¥=V¥2,m)=e/r, 7° = (R@)’ + & — a)’, (3.18) 
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l ° 
' ,1ld 
M(z) = uke | ER*(1 + ay" . a dz, 
ly 


R’ = dR/dz. 


Since confusion with regard to partial differentiation has now been removed, we have 
restored the usual notation, y being regarded as a function of z and z . It should be 
noted that on account of the occurrence of R(z) in r, the kernel K(z, z)) is not symmetric. 
Since # = 0 outside the gap, the limits of integration in the expression for M(z.) may 
be replaced by those corresponding to the ends of the gap. 

4. Some radiation problems. In setting up simplified mathematical models for the 
treatment of radiation problems, it is essential that we carry into the model sufficient 
conditions to yield a definite mathematical problem, without putting in so many condi- 
tions as to make the problem insoluble. The following discussion is based on the fact 
that a field of outward radiation is determined by the assignment of E, over a surface 
S which forms the inner boundary of the infinite region. 

The mast antenna. The common problem of macrowave radiation involves a vertical 
mast above a horizontal earth (Fig. 3). Let us make the usual assumption of infinite 
conductivity. We think, then, of an infinitely conducting cylinder AB above an in- 
finitely conducting plane. The apparatus is fed by a generator D, the terminals of which 


A 


- 
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Fia. 3. The mast antenna 
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are connected to AB and to the plane at C. We have then the problem of determining 
a Maxwellian field satisfying the following conditions: 
(i) outward radiation at infinity, 


(ii) E, = 0 over AB, over the plane earth, and over the connecting leads. 
Where does the energy come from, which reaches infinity? It cannot come from 
the perfectly conducting surfaces, since E, = 0 implies that the Poynting vector is 


tangent to a perfectly conducting surface. Hence the only possible source of the energy 
is a surface drawn around the generator D. (If that surface were perfectly conducting, 
the generator would be short-circuited, and there would be no radiation at all.) The 
rest of the apparatus serves only as a guide to the radiant energy. 

It would appear therefore that we can obtain a determinate problem in the case 
of Fig. 3 only by assigning E, over a surface which bounds the generator. If this were 
done, we would have E, all over a surface consisting partly of the aforesaid surface 
and partly of the surfaces of the conductors, and so we would have a determinate 
problem. 

But in the actual treatments of such problems, no reference is made to the geometry 
of a surface bounding the generator, nor even (usually) to the geometry of the leads. 
A discussion along the lines suggested above would be, apparently, unnecessarily com- 
plicated from a practical viewpoint. It is usual to simplify the problem by substituting 
a model as shown in Fig. 4 for the actual system. Now the generator D is reduced to 


oe 





Lz D 








Nel 





Cc 
Fic. 4. The mast antenna with gap. 
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small proportions and placed between the antenna and the ground. To provide a surface 
for boundary conditions, the geometrical surface of the antenna is continued down past 
the generator to the ground. 

It is not, of course, suggested that the above arrangement is ever employed in 
practice. But it is the arrangement employed in mathematical discussions, whether 
the authors admit it or not. The hope is that the behavior of the model shown in Fig. 4 
(deduced mathematically) will agree with the behavior of the actual apparatus of Fig. 3. 
On this point we have no assurance from theory. 

Figure 4 introduces the gap; by this we mean the surface BC, bounding the generator. 
Many writers have proceeded as if antenna problems could be discussed without reference 
to a gap. To do so, however, is to leave out the essential part of a radiating system. The 
gap is the only source of radiant energy. 

Now if we knew E, over the gap, we would have a determinate problem in Fig. 4. 
But we do not know E, over the gap. This may be expressed by saying that we do not 
know the (electrical) structure of the gap. But if the gap is short, its structure does not 
seem to be important. It is then possible to disregard the structure, and express the 
radiation properties of the antenna in terms of only one characteristic of the gap, namely, 
the voltage or potential difference. This is the integral of E, taken along BC, with suitable 
sign. 

The apparatus used by Brown and Woodward’ for the experimental measurement of 
impedance is similar to that shown in Fig. 4, except that their gap is shorter and the 
mast is fed by a coaxial line which comes in from below through the ground. 

The cigar antenna. It is customary and legitimate to eliminate the earthplane in 
Fig. 4 by a method of symmetry. This leads to the cigar antenna (Fig. 5). This self- 


G Ss 


C ZL » 


Fig. 5. The cigar antenna. 














contained radiator is a mathematical fiction, but it could be realized in practice. We 
think of a hollow copper cigar, cut in two to leave a gap. Inside there is a generator, 
with terminals connected to the two ends of the cigar. (The two generators obtained 
by reflection of Fig. 4 in the earth are replaced by one generator in Fig. 5; but this is a 
triviality, because we are not interested in the interior of the cigar. What interests 
us is E, over the gap G.) 

Although the source of the cigar antenna is to be found in the application of the 
reflection process to the mast antenna, it suggests a generalization in which symmetry 
is abandoned (Fig. 6). Figure 7 shows a coaxial line terminated in a spherical knob. 
The model for this (Fig. 8) is the limiting case of an asymmetric cigar, one end being 
reduced to a small circular patch on a conducting sphere, surrounded by a ring gap. 

The later detailed applications in this paper deal with a cigar antenna, in general 
asymmetric. 

The microwave antenna terminating a coaxial line. Let us now consider a microwave 


8G. H. Brown’‘and O. M. Woodward, Proc. I. R. E., 33, 257-262 (1945). 
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Fic. 6. The cigar antenna (unsymmetric). 
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Fig. 7. Coaxial line terminated by a spherical knob. Fig. 8. Spherical antenna with ring gap. 
apparatus. Figure 9 shows a coaxial line, terminated in an antenna. We may introduce 
the idealization that the line extends to infinity. 

Let us see whether this apparatus gives a definite problem. For surface S we take 
the surface of the antenna and the core of the coaxial, together with the inner and outer 
surfaces of the sheath of the coaxial. Over this surface we have E, = 0. Conditions at 
infinity, however, are not obvious. We cannot say that there is outward radiation in all 
directions—there is nowhere for this radiation to come from. We expect that there will 
be outward radiation everywhere except in the neighbourhood of the coaxial line, and 
that inside the line there will be inward radiation (together with some outward radiation 
if there is incomplete tuning). 

Thus, while Fig. 9 presents difficulties as far as the definition of a determinate mathe- 
matical problem is concerned, these difficulties are not those encountered previously. 
There is no troublesome gap—unless we like to say that the gap is at infinity. 
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Obviously some simplification is necessary for mathematical attack. The usual 
simplification is drastic. We retain only the horizontal-arms of the antenna, wiping 
out all the rest of the apparatus. We join the arms (dotted indication in Fig. 9) and 
regard this cylindrical boundary as a gap. In fact, we return to the cigar antenna! It 
is indeed remarkable that impedance calculations obtained from such a simplified model 
should have the practical value that they do seem to possess. 

The idealization of the microwave antenna of Fig. 9 to the cigar antenna of Fig. 5 
involves us in the old difficulties concerning the gap. This seems particularly objection- 
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Fic. 9. Two-arm antenna terminating coaxial line. 


able because in actual antennae the gap is not very small compared to the length of 
the antenna. Is there not some way to discuss the problem without introducing a gap 
artificially? It appears difficult in the case of the arrangement of Fig. 9. 

With a different form of apparatus, however, the situation is more hopeful. Figure 10 
shows a coaxial line, with the core projecting as an antenna, while the sheath spreads 
out into a conducting plane. The surface S consists of the surfaces of the core and 
antenna, with the inner surface of the sheath and the surface of the conducting plane. 
Over S we have E, = 0. On the infinite hemisphere above the antenna we assign out- 
ward radiation. At infinity down the line we assign a direct wave (amplitude C,) and 
a reflected wave (amplitude C,). The problem of determining the ratio C,/C, is then a 
definite one. We approach the important problem of reflection without using the artificial 
concept of impedance. 

To bring the preceding problem more into line with what was said earlier regarding 
determination, we may change the conditions down the line. Let us say that E, is assigned 
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on a cross-section a long way down. Now we have a definite problem, since E, is given 
over the bounding surface, with a condition of outward radiation in the hemisphere. If 
we find H, on the section of the line, we can find the ratio of the reflected to the incident 
wave. 

The case of Fig. 10 should therefore be treated with considerable respect, because 
it represents an actual situation in which the radiation problem is definite; it is un- 
necessary to simplify the apparatus drastically for mathematical purposes.” 
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Fic. 10. Coaxial line with antenna projecting 
beyond conducting plate. 


Before leaving the case of Fig. 10 we may mention how this apparatus may be con- 
sidered as of the gap type. If we consider only phenomena above the plane of the flat 
conducting surface, we have a situation very like that of the conducting mast; the only 
difference is that the gap is now around the base of the mast instead of being below it. 
Once more we may appeal to the method of symmetry and discuss the problem as a 
cigar problem (Fig. 11). Now the surface S consists of a complete cigar with a. disc 
attached at its center. We must assume E, given over this disc, which is in fact the gap 
in a new form. We are not to expect that the field has continuous derivatives if we pass 
through the disc. We may call Fig. 5 a cigar-with-a-band and Fig. 11 a cigar-with-a-frill. 
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Fig. 11. Cigar with a frill. 


*Experimental results for this case have been given by Brown and Woodward, loc. cit. 
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The antenna in the wave-guide. Figure 12 shows an antenna in a wave-guide. This 
problem lends itself to the same satisfactory mathematical formulation as for Fig. 10. 
There are only two differences. Instead of a hemisphere at infinity above, we have the 
opening of the wave guide at infinity on the right, and there we impose a condition of 
outward radiation. (That is in the case of transmission; for reception, we must allow 
for a direct and reflected wave, with only a wave down the coaxial line.) Secondly, the 
five faces of the guide (if it is of rectangular section) increase the complexity of the 


geometry. ° 
Should we wish to introduce the gap concept, we may use the method of images 
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Fig. 12. Antenna in a wave-guide. 


when the guide is rectangular. This means that we forget all about the walls of the 
guide and the coaxial line, and consider a doubly infinite pattern of cigars-with-frills 
in empty space. 

These problems of radiation have been discussed in general terms at some length, 
because such discussions have unfortunately been lacking in the literature of radiation. 
It seems essential that we should have at least an idea of the conditions that make a 
problem definite from a mathematical viewpoint, and also that we should recognize the 
simplifications that are made to reduce practical problems to a form that admits mathe- 


matical treatment. 
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THE GENERAL PROBLEM OF ANTENNA RADIATION AND THE 
FUNDAMENTAL INTEGRAL EQUATION, WITH APPLICATION TO 
AN ANTENNA OF REVOLUTION—PART II* 
BY 
J. L. SYNGE 
Carnegie Institute of Technology 
5. The gap. Let us assume the gap to be cylindrical, of radius R = a; let it extend 
from z = —7toz = 7. Let us suppose that on the gap 


E = —3V/n, (V = const.) (5.1) 
so that V is the potential difference across the gap. Then, by (3.18), 


"1d a (5.2) 


M(&) = mck*VN(q), Ne) = — rd 


Now 


—1l ékr 


r=a’+(z—xZH)’, y=r'e’, 


ldy_§ 1 Fiijs , 
Me he 9, + Oa + ix), 





r® 
x = ie (1 — 3k’r? — cos kr) + “ns (kr — sin kr), 6.3) 
k®r® k rT 
% = a (kr — sin kr) — in (1 — cos kr). 
a oY k’r 
Note that x; , x2 are power series in positive powers of kr. We obtain at once 
. 1 n — Zo n + 2 
N (Zo) = ZY ~ 271/ + 9 } 
kn \[a® + (n — 20)*)'? [a + (9 + &)")"” 
1 k’a’ ' 2 241/2 
+ 5 7 {In k(n — 20 + [a + (n — &)]”) 
(5.4) 
+ In k(n + 2 + [a? + (n + 2)*)'”) — In k’a"} 
— = Le + tx) de, 
We shall now make two important simplifying assumptions. The first is 
ka K 1. (5.5) 








*Received March 31, 1947. This paper is a continuation of the paper by G. E. Albert and J. L. Synge 
which appeared in this Quarterly 6, 117-131 (1948). 
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This means that the radius of the antenna at the gap is small compared with the wave 
length. Then we have approximately 


1 ” — Zo n + 2 \ 
N y =, { =e 5.6 
© - in iar + QT? | + eT saa 


provided a/7n is not large. The second assumption is 


a/n <1. (5.7) 








This means that the gap is long compared with the radius of the antenna at the gap. 
Then (5.6) gives approximately 


_— 7 n — %o ate} 
Me) = pt + EST, “ 





and so 


NG, = é for | | < 9, 


(5.9) 
N (Zo) = 0 for | Zo | > 7. 
Substitution in (5.2) gives 
M(z.) = 2axckV/n for | Z| < 7 (in gap) 
(5.10) 
M(z) = 0 for |2| > (outside gap) 


6. Impedance and relative current. In this section we introduce the impedance Z 
and the relative current ¢(z), and show how Z is found when ¢(z) is known. The argu- 
ment is exact; we understand by N the exact expression (5.2) rather than the approxi- 
mation (5.9). 

We define the impedance of the antenna to be 


Z = V/I(0) (6.1) 


It may seem unnatural to use (0) in defining impedance. The point z = 0 is at the center 
of the gap, and there is no current there. In fact, /(0) means 2rcaH (0), according to 
(3.14). It might seem better to use the currents at the ends of the gap. But it appears 
simpler to use J(0) as basic; we can easily pass to the other definitions, if required. In 
the case of a very short gap (Sec. 8), these subtle distinctions disappear, for we have 
then 
I(—n) = 1(0) = I(n) 
approximately. 
Let us write 
g(z) = I(z)/I(0); (6.2) 


this will be called the relative current. 
On dividing (3.17) by I(0), we get 


ls 
Kz, 20)6(z) dz = imck’Z N(z), d<a<& (6.3) 
i; 
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With this equation we associate the boundary conditions 
¢(0) = 1, ¢(l,) = o(L) = 0. (6.4) 


The equation (6.3) contains the unknown function ¢(z) and the unknown constant 
Z. If we knew ¢(z), we could calculate Z at once, giving any value to 2, . If we have only 
a rough idea of ¢(z), it would be better not to take a definite value of z) , but to introduce 

a. weighting factor f(z)), and calculate Z from 
Z = — fi? Leo) deo i? K(z, 20)(2) dz (6.5) 

- Si? N@o)f Go) dzo 

For the present let us leave the weighting function f(z) arbitrary except for the 
assumptions that it is continuous, has a continuous derivative, and satisfies the end 





conditions 


fil) = fil) = 0. (6.6) 


Let us write J for the numerator in (6.5) and understand the limits of integration 
(l, , 2). Then 





J = [ff Ke, «o@se) de de 
(6.7) 
a . 
= — [fe cesta) de des +H ff We, ed @S 60) de den, 
and on integration by parts 
I= [ff ve, zdx(e, 20) de des 6.8) 
where 
x(z, 20) = —o'(z)f"(20) + k’o(z)f (20). (6.9) 
Now we may write (6.8) as follows: 
J=JS,+Jd2, 
l 
J; = Parties 2) dz dzo ; 
(6.10) 





J. = I/ ~ {x(z, 20) exp (tkr(z, 20.)) — x(z, z)} dz dz , 


r(Z, Zo) 
r 2 2 2 
[r(z, 20)] [R(z)}’ + (@ — 20)’. 
For a thin antenna, r is small for z = 2) . However, the integrand in J; remains finite 
as R — 0. This is the reason for splitting J as above, and forms the basis of later ap- 


proximations. For the present the argument remains exact. 
The integral J, gives 


Il 


J, = Ju t+ Jia + Ji, (6.11) 
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where 
Is 
Jy 2 L x(z, 2) dz, L = —In(ka’), 
ls 
le a’ ° 
Tae / “xe, 2) In ow de (6.12) 


Is 
Jis = / x(2, 2) In {H’[l, — 2 + (R’ + (LL — 2)")'"Ile — Le + (R’? + & — 1))'")} az. 


Note, for later approximation, that J, , J; remain finite for an infinitely thin antenna. 
We have, by (6.6), 


Ju =L | {-#'@s'@ + Pees} de 


#1, 
(6.13) 
als 
=L] {f"@) + PI/@\6@) de. 
Now (6.5) may be written 
a. Jil, f) + Ji, f) + Jis@, f) + Jo, f) (6.14) 





ack” i? N(z) f(z) dz 


This notation puts in evidence the dependence of the J’s on the two functions—the 
relative current ¢ and the weighting function f. 

The function f is at our disposal. We see from (6.13) that J,, would vanish for an f 
sinusoidal in kz. However, unless the whole length of the antenna is a multiple of 3), 
there exists no such function with continuous derivative, satisfying the end conditions 
(6.6). So we approach a sinusoidal f by a limiting process, in which (on attaining the 
limit) the continuity of the derivative is lost, but (6.14) remains true. 

Let ¢« be any small positive number. We define a function f,(z, €) as follows: 


lL<z2< —e: K(f, -) = —sin kl, sin k(l, — 2), (6.15a) 


—e<z<e :K(6fil(z, 6) = Ke) cos kz + §sin kzsin k(l, + 12) 
(6.15b) 


— 4 cosec ke sin k(l, — 1,)(1 — cos ka), 


II 


e<xz2< :KO)f:C@, © —sin kl, sin k(l. — 2), (6.15¢) 
K(e) = —sin kl, sin kl, — 4 tan $ke-sin k(l, — 1). (6.16) 
To avoid complicating the argument, we assume 


ain kl, = 0, sin kl, ~ 0. (6.17) 


This means that neither arm of the antenna, measured from the center of the gap, is a 
multiple of 4\. Such critical cases must be approached by a special limiting process. 
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We note that f,(z, €) is continuous, with continuous first derivative, and satisfies 


z|>e:ffz,) + Ffilz, 2 = 0, (6.18a, ¢) 


_F sin k(l, — 1) 








‘gi : f;’G, €) a k’ f(z, €) = 2K(e) ie the (6.18b) 
Then, by (6.13), in an obvious notation, 
~ k? sin k(l, sensi L,) . 
J33:(¢, file) = —- 2K(e)sin ke _ o@) dz. (6.19) 
Let us write 
f(z) = lim file, €), (6.20) 
so that 
= s< 2: fi(z) = sin k(l, — z)/sin kl, , (6.21a) 
OS <2 = ib: fi(z) = sin k(l, — z)/sin kl, . (6.21¢) 


Proceeding to the limit « — 0 in (6.19), we get (since ¢(0) = 1) 


sin k(l, — 1) 
Jno, f;) =kLT, T= 
(dy J ; sin kl, sin kl, (6.22) 
Let us now put f,(z, €) for f(z) in (6.14) and proceed to the limit « — 0. In this limiting 
process, f,(z, €) and its first derivative remain finite in the whole range (I, , 12), and so 
the contributions to J,., J:, , J2 from the range (—e, €) vanish in the limit. Thus we 





get 
} 4 =_ — 5 kLT + J ,(¢, L) = J13(¢, ti) + J (¢, fi) (6.23) 
eh fit N@f@) de 
where f,(z) is as in (6.21) and I as in (6.22); L = —In(k’a’). 


This is accurate, and indeed holds for the general N of (5.4) as well as the more 
particular N of (5.9). The fact that f,(z) has a discontinuous first derivative at z = 0 
creates no trouble. It is of course understood that in evaluating J, by (6.10) and J, , 
Ji, by (6.12), we are to put 


x(Z, Zo) = —@'(z)f1(zo) + k’b(2)f; (20). (6.24) 
On substituting for N from (5.9) in the denominator of (6.23), we obtain 


t 





Z = 7-75 (RLY + Jusld, fr) + Sis, fr) + Je, fd}, 
(6.25) 
_ Sinkn , ,,1— cos kn 
"key ; kn 


This is the formula we shall use in the later work. It is accurate except for the approxi- 
mation involved in (5.9). ; 

Note: Do not confuse the range (—e, €) with the gap (—7, 7). The former is merely 
a mathematical device, introduced to eliminate ¢ from the first term in the numerator 
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of (6.23). As a matter of fact, we shall make no further use of this ¢; it has done its 
work in providing the formulae (6.23), (6.25). 

7. The thin antenna. We cannot expect to get results for an antenna of general 
form without a considerable amount of calculation. But if the antenna is thin, the large- 
ness of ZL may be used as a basis of approximation. In the present section we obtain 
the principal parts of the current and impedance for a thin antenna. This current is, 


VA i 


7 yf. 








— 

















Fia. 13. Cylindrical antenna with spheroidal ends. 


in fact, the familiar sinusoidal distribution. The present derivation of this distribution 
may be of interest because previous derivations have been by no means clear. More- 
over, in the present method, it is not necessary to assume-that the antenna is cylindrical; 
it is merely necessary that the radius R be small throughout. The ends of the antenna 
require no special treatment. Further, our formula for impedance contains a shape 
term. The method of the present section does not open up a process of successive ap- 
proximations; that will be given in Sec. 10. 

It will be well to mention here assumptions which will be introduced explicitly later: 


The antenna is thin at the gap (ka < 1). (7.1) 
The antenna is thin throughout (kR(z) « 1,1, < z < l,). (7.2) 
The gap is long compared with the radius (a/n < 1). (7.3) 
The gap is short compared with the wave length (kn < 1). (7.4) 


Obviously (7.2) contains (7.1). 
We turn back to the exact formula (6.14), in which f(z) is arbitrary. By (6.13) we 
have 


Iu =L | (6"@ + KONO ce. (7.5) 


1 


We rearrange (6.14) in the form 
la 
[ (o"@ + K6@) — inc? ZL*N@)$@ de = -—L"VJut Js + Js). (7.6) 


Now make the assumption (7.2). Then L is large and the right hand side of (7.6) is 
small of order L~*, provided ¢(z) and ¢’(z) remain bounded as L tends to infinity. Since 
f is arbitrary except for the end conditions (6.6), it follows from (7.6) that 


$'"(2) + K’o(2) — imck*ZL"N(2) = L“5(2), (7.7) 
where §(z) is some finite function. Integration gives 


o(z) = a cos kz + Bsin kz + imckZL"' N(t) sin ke — #) dt + L™'S(z), (7.8) 


where G(z) is some finite function. This solution is subject to the three conditions (6.4), 
and if we knew G(z) we could find a, 8, Z. But we do not know G(z), and can merely 
make use of the fact that the last term is of order L™*. 
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Whatever G(z) may be, the three equations for a, 8, Z are consistent provided 


qa 
sin Bl, | N(t) sin k(l, — #) dt 
0 
x 0. (7.8a) 
la 
sin kl, N(4) sin k(lp — t) dt 


Let us make the assumption (7.3), so that we may use (5.9) for N. Then the above 
condition for consistency reads 


(1 — cos kn) sin k(l, — 1) + 2 sin kn sin kl, sin kl, ¥ 0. (7.8b) 


Assuming that kn, kl, , kl, are such that this inequality is satisfied, we obtain a, B, Z 
from (6.4), and substitution in (7.8) gives accurately 


$(z) = i(z) + L'Q(2), (7.9) 





$,(2) 
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Fic. 14a. Infinitesimal gap at center. 


where Q(z) is an unknown finite function and ¢,(z) is given by 


l, <2< —n : K(n)¢,(z) = —sin Kl, sin k(l, — 2), (7.10a) 
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Fic. 14b. Finite gap (n = 1/4) at center. 


—n <2 <7 : K(n)bi(z) = K(n) cos kz + 3 sin k(l, + 1.) sin kz 


(7.10b) 

— 4 cosec ky sin k(l, — 1,)(1 — cos kz), 
n<2z<l,. :K(n)¢,(z) = —sin kl, sin k(l, — 2), (7.10c) 
K(n) = —sin kl, sin kl, — } tan $k sin k(l, — |). (7.11) 


(These formulae should be compared with (6.15), (6.16). Note that /, — 1, is the length 
of the antenna.) 

Physically, ¢,(z) represents the principal part of the relative current for a thin 
antenna with a gap which is long compared with the radius of the antenna, and becomes 
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a better approximation the thinner the antenna becomes. Graphs of ¢,(z) are given 
in Figs. 14a-g (for discussion, see Appendix, p. 155). Outside the gap, ¢,(z) has the 
sinusoidal distribution, so basic in antenna theory. 

To get the impedance, we substitute for ¢ from (7.9) in (6.25): This gives 


* ~o {kKLT + J i2(¢1 » fi) + isi ’ fi) + I$: , fy) 


(7.12) 
+ L"[Ji2(Q, fi) + Jis(2, f.) + J2(Q, f.)]}. 
Here everything is known, except ©. 


Let us sum up: 
For a thin antenna with a gap much longer than its radius, the relative current is given 


by (7.10), with an error of order L~*, and the impedance is 


Z; - 7s {kLT + J i2(or ’ fi) + Ji3(¢: » fi) + J (1 ’ fd}, (7.13) 
TC. 





(2) 
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Fia. 14c. Finite gap (7 = 1/2) at center. 


ale 


. -1 
with an error of order L’’. Here 





L = —In(k’a’), a = radius at gap, k = 2x/d, 
_ sinkyn , ,,1— cosky a : 
h =— kn + oT kn ’ ( n < 2 < n) 1S gap, (7.14) 


sin k(1, — L,) —_ ae % 
sie tell ala Lana” cot kil, cot ki, . 


r= 
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We must remember that f, is given by (6.21) and does not depend on e, whereas ¢, is 
given by (7.10) and does depend on 7. 

On account of the assumed thinness of the antenna, we can simplify the expressions 
(6.10) for J, and (6.12) for J,; . In fact, we shall put R = 0 in J, and J; , since the 
consequent error is of order ka, and so is negligible in comparison with L~*. Thus we 


write 
als a’ 
Ju=] x2) ng ae, 
ta P 
Ji3 = [ x(z, 2) In [4k’(l. — z)(z — 1,)] dz, 


— (7.15) 


ls ls 
J, = [ [ | z — 20 |" {x(@, 20) exp tk | z — 2 | — x(@, z)} dz dz, 
ej i 


x Se 


—$i(z)fi(eo) + k’oi(2)f (20). 


x(z, Zo) 


We note that the shape of the antenna is now involved only in J,, , and that there is 
no contribution to J,. from cylindrical parts of the antenna. 
Putting Z, = R, — iX, , the approximate resistance and reactance are (with an 


error of order L~’) 





1 
R, _ Amckh J 22 ’ 
(7.16) 
1 . 
xX, = javkh (ALT + Jiz + Jis + J), 
where 
Ja + tdeg = Je. (7.17) 


Thus the resistance is independent of shape. 

8. Thin antenna with short gap. The calculation of the impedance from (7.13) is 
direct. The result will, of course, depend on the gap-length 27, and will be very involved 
on account of the complexity of (7.10). Let us therefore, for simplicity, introduce the 
assumption (7.4); we consider the gap short compared with the wave-length, but still 
long compared with the radius of the antenna. 


If kn is small, h = 1 approximately. Further, by (7.10), we have approximately 
lL<2z<0: ¢;(z) = sin k(l, — z)/sin kl, , (8.1a) 
O<2<b: ¢,(z) = sin k(l, — z)/sin kl, . (8.1¢) 


(Note: the consistency condition (7.8b) becomes sin kl, sin kl, # 0 for small k7.) 
On comperison with (6.21), we see that ¢, = f, . This greatly simplifies the work by 
introducing symmetry. We have 


x(Z, Zo) = —i(z)b! (Zo) + k*d,(z)br (20), (8.2) 
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Fig. 14d. Infinitesimal gap at point dividing antenna in ratio 3:5. 
and so 
lL <2<0:x(,z) = —k’ cosec’ kl, cos 2k(l, — 2), (8.3a) 
0<2z< lI, :x(z,z) = —k’ cosec’ kl, cos 2k(l, — 2), (8.3) 
z< 0,2 <0:x(z, 2%) = —k’ cosec’ kl, cos k(2l, — z — z), (8.4aa) 
z2>0,z% <0,orz<0,%>0: 
(8.4ac) 
x(z, 2.) = —k’ cosec kl, cosec kl, cos k(l, + l, — z — 2), 
z> 0,2 > 0:x(z, 2%) = —k’ cosec’ kl, cos k(2l, — z — %). (8.4¢cc) 
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It is best to combine J;3; and J, in (7.15). We have 


ll |z — 20 |~*x(z, 2) dz dey 
(e) 


i. | " xe, 2) Inf4k*(y — 2)(e — 1,)] de — Bin(2ke) : adele 


a 


(8.5) 


Here f {,.) means integration over the square (J, , 1.) with omission of the strip | z — z | < 
e«. Now 


J, = lim I |z — 20 |” '{x(z, 20) exp tk | z — 2 | — x(z, 2)} dzedz, (8.6) 
(e) 


«0 


and so 


Ji3 + J2 = lim ff |z — 2 |~'x(z, 20) exp tk | z — 2 | dz dz, 
(e) 


«-0 


(8.7) 
+ 2in(2ke) | x(z, 2) de |. 





$, (2) 
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l4e. Finite gap (n = 1/4) with center dividing antenna in ratio 3:5. 
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This double integral is best evaluated by rotating the z, 2, axes through half a right 
angle in the plane of the integral. It is convenient to introduce the function 

In y = 0.5772. (8.8) 


P(x) = | c—4 dt = Ciz+iSiz— Inyz, 


We find 
(l, — 1)? 


¢) ln Ss 


(Ji3 + J2)/k _ (C2 = 4k? 7? 


+ (4k) {c, a 4(€:C2 + 1)} 
(8.9) 


— B(2kl,) {co — c, + te2(co — ¢;)} 
—_ ©(—2kl,) {c. i 1; (C2 _ ¢:)}, 


Cc. = cot kl, , c, = cot ki, , 21 = 1, — l, = length of antenna. 


(2) 
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Fic. 14f. Infinitesimal gap at point dividing antenna in the ratio 1:3. 


When we substitute from (8.9) in (7.13), we get the following expression for the 
principal part of the impedance of a thin antenna, expressed in ohms (1 Heaviside unit = 


120 zc ohms): 
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Zi = R, = 1X, 
sin 2kl war: 
= 80 FTL em EL {2 In | “2 | — @(4kI)(i — cot 2k) 
(8.10) 
+ &(2k | 1, |) — cot k | i, |) 
+ &(2k/,)(i — cot uy} — iX,, 
where 
, 0 a 
X, = 30J,/k = —30k cosec’? kl, J ~ ¢08 2k(, — 2) In ape de 
(8.11) 


2 
a 


ls 
¢ cosec’ k 30s 2k(l, — are (1, 
30k cosec” kl, J cos 2k(l, — z) In (Ro) dz 
Here z = |, , 2 = I, are the ends of the antenna (1, < 0 < l,); the gap is short and at 
z = 0; 2/ is the length of the antenna; a is the radius at the gap, and R(z) the radius 
at the general point; k = 22/i, where is the free wave length. 
If we take the gap at the middle, then —l, = 1, = l, and (8.10) becomes 


Z, = 60 cot kl{2i In(l/a) — (4kl)(¢ — cot 2kl) 
(8.12) 
+ 26(2kl)(i — cot kl)} — iX,. 


Except for the shape term X, , this agrees with the formula given by Brown and King’ 
using the method of Labus.’ There is also agreement with the principal part of Schel- 
kunoff’s formula’ and with the formula of Hallén.* Schelkunoff’s treatment: of the 
influence of shape is difficult to follow. Hallén includes a shape term in his equation 
(26), but later specializes to a cylindrical antenna, so that a formula such as our (8.11) 
does not occur explicitly in his equation (39). Owing to the inadequate treatment of 
the gap in the work of Hallén and Schelkunoff, the validity of their higher approxima- 
tions is open to question. It must be remembered that an error of the order L~* is ad- 
mitted in our formulae (8.10), (8.12). ; 

In (8.12) it is not assumed that the antenna has z = 0 for equatorial plane of sym- 
metry. Deviation from this symmetry influences*X, , but not the other terms. 

Let us consider an antenna with the gap at the center and total length nearly 3), 
so that 


kl = jr +, (8.13) 
where ¢ is small. Then, approximately, with errors of orders e° and e’ respectively, 
cot kl = —e, cot kl cot 2kl = —}. (8.14) 


1G. H. Brown and R. King, Proc. I. R. E. 22, 457-480 (1934). 

2J. Labus, Hochfrequenztechnik und Elektroakustik 41, 17-23 (1933). 
38. A. Schelkunoff, Proc. I. R. E. 29, 493-521 (1941). 

‘E. Hallén, Nova Acta Reg. Soc. Sci. Upsaliensis 11, No. 4 (1939). 
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Then we may neglect parts of (8.12) and write 
Z, = R, — 1X, = —120%e n* — 306(27r) — iX,, 


R, = 30(log 2ry — Ci 2x) = 73.138, 
(8.15) 


X; 


120¢ In + X, + 30Si 27, 


1206¢ In © + X, + 42.54. 





$ (2) 
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Fic. 14g. Finite gap (» = 1/4) with center dividing antenna in the ratio 1:3. 
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The case where the antenna is approximately of total length 3, but with the gap 
not at the center, is also of interest. Again we have (8.14), and substitution in (8.10) 
gives approximately 





Z, = —30 cosec’” kf 62) + 4ie In | ie } — 1X,, 
R, = 30 cosec’ kl,(In 2ry — Ci 2r) = 73.13 cosec’ kl, , ' (8.16) 
X, = 120¢ cosec’ kil, In i | + X, + 42.54 cosec’ kl, . 





Thus, by moving the gap away from the center of the antenna, we increase the resistance, 
but the tuning of the antenna to make X, = 0 is more difficult, because the derivative 


of X, with respect to ¢ is greater. 

We see that the problem of matching the antenna to a coaxial line, as far as reactance 
is concerned, depends on the shape term X, in an important way. This term is discussed 
in the next section. 

9. The shape term in the reactance. Let us consider the term X, , given in (8.11). 
For simplicity, let us assume that the gap is at the center of the antenna and that the 
antenna has an equatorial plane of symmetry (z = 0). Then (8.11) reads 


l 2 
X, = —60k cosec’ kl / cos 2k(l — z) In Rey dz, (9.1) 


2] = length of antenna. 


We note that X, receives no contribution from cylindrical portions of the antenna. 

We can now settle the vexed question of contribution from the ends of a cylindrical 
antenna, by supposing the antenna to be a cylinder (R = a), terminated by spheroids 
of semi-axis b (Fig. 13, p. 138). Then on the ends 





(R@) , @—-l+bdP _ 
7+ 7 = 1. (9.2) 


No assumption is made at first about the magnitude of b. Equation (9.1) gives 
: ; . 
X, =-60k cosec? kl J cos 2k(l — z) Inf{1 — (¢ — 1 + b)?/b?] de. (9.3) 
i-b 


Since the logarithm breaks into the sum of two logarithms, this integral is easily evaluated 
in terms of Ci and Si functions; we find 


X, = 30 cosec? kl{sin 4kb(Ci 4kb — Ci 2kb) 
(9.4) 


— cos 4kb Si 4kb — (1 — cos 4kb) Si 2kb}. 
If kb is small, this approximates to 


X, = 120 kb. cosec’ kl. (In 2 — 1) ohms. > (9.5) 


Since this tends to zero with kb, we see that there is no contribution to reactance from 
the ends of a cylindrical antenna cut off square at the ends. (Of course, “contribution from 
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the ends” implies some mathematical division of reactance into “contributions” of 
various sorts; our statement refers to the division we have made in using the integral 
(9.1).) In fact, X, = 0 to our order of approximation if the ends are rounded for a length 
b comparable to the radius a of the cylindrical part even for non-central gap. 

Spheroidal antenna. The impedance of a spheroidal antenna can be calculated ac- 
curately by means of spheroidal functions.” However for a thin spheroid, the present 
method may be used. Let us take the gap at the center. Then 


(R(z))*/a® = (P — 2°) /P, (9.6) 


and so by (9.1) 


l 
X, = 60k cosec® kl | cos 2k(U — 2) Inf(P — 2*)/P] de 
0 
= 30 cose? ki{sin 4k1(Ci 4kl — Ci 2kl) (9.7) 
~ cos 4kU(Si 4kl — Si 2k) — Si Qkl}. 
For 1 = /4, kl = 42, we have 
X, = —30 Si 2r. (9.8) 


Referring to (8.15) with « = 0, we see that X, = 0; the reactance of a thin spheroidal 
half-wave antenna (with the gap in middle) is zero. This fact is mentioned by Schelkunoff 
(loc. cit.), but the reason for this statement is not clear. 

Conical antenna: For a symmetrical thin conical antenna fed at the vertex, we put 


R(z) = Bz (9.9) 
where 6 is the semi-angle of the cone. Equation (9.1) gives at once for the shape term 


. X, = 120 cot kl In £ + 60. (9.10) 


For the approximation to be valid, we should take 6 of the same order as ka. 

Graphs of impedance are given in Fig. 15 (see Appendix). 

10. Successive approximations. The method of Sections 7 and 8 gives the principal 
part of the impedance for a thin antenna, but it does not open up a method of successive 
approximations. To get such a method, let us return to the exact integral equation 
(6.3) and write it in a slightly different notation as follows: 


pls a2 ad le 
~j* a do at +i? [” v(t, 26 dt = inck’ZN@,  (h <2 <1). (10.1) 
Ji, OZ © 2s 


Let us introduce the integration operator S such that 
Si@ = | fo at. (10.2) 
/o 


5J. A. Stratton and L. J. Chu, J. Applied Physics 12, 241-248 (1941); L. Infeld, Q. Appl. Math. 5, 
113-132 (1947). 
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15a. Resistance of thin antenna with infinitesimal gap at point of quadrisection. 
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on account of (6.4). Here A is a constant, 
Is 
A= | v(t, 00" dt. (10.4) 
ty 
Now we operate on (10.1) with S’; this gives 


S| wt," dt — Az + BS? | " W(t, 2)o(t) dt = ixck®ZS*N(e). (10.5) 
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Fig. 15b. Reactance of a thin cylindrical antenna with infinitesimal gap at point of quadrisection. 
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To treat these integrals, consider the operator B, defined by 
la 


Bf(z) = S v(t, z)f(t) dt 


I] 


als 


[soul vena 


II 


ll 


B, f(z) + B.f(2), 











where 
p* dy —_ 2 2 
Byje)= | sae} 2, it, WP = (ROP + W - 0, 
Pe he ** exp ikr(t, y) — 1 
B.f(z) = | f() dt i ge dy. 
Chen 
ase ~ [wetter te ol nga 
Fs, —t+ {(R()) +.°}'7 
2 als 
= (| + | in k[z — t+ {(R(D]? + @ — 07}'7)F(O at 
~ (/ +] ) In k[—t + {(R(OY + 2371/0) at 
“1, “0 
z pls 
= (/ _ | ) im Aij2—t| + {(RO? + & — O77 (0 dt 
re (| - [”) In kt | ¢| + {(ROP + Cy?) dt 
ls nla 
+ [ nfer(oy’s® at — | nfaR(ors(o at 
= B.f(z) + B.f(2) — In(k’a")Sf@), 
where 


B.f(z) = ([ - |) Inkl |z2—2t| + {{[R(@)? + & — dy? 1f( dt 
_ (| _ | ) Inkl | ¢|) + {(R@? + PY 7 VO dt, 


Bj) = | Inla/R(O]*f() at. 


al | 
In the notation of (10.6), (10.5) reads 


Bo'(z) — Az + k°SBo(z) = ixck’ZS°N(z2), 


(10.6) 


(10.7) 


(10.8) 


(10.9) 


(10.10) 
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or, since S¢’(z) = ¢(z) — 1, 


L{o(z) — 1 + k’S’o(z)] + Fo(z) — Az = trck’ZS°N(z), (10.11) 
where L = — In k’a’, and F is the operator 
F = (B,+ B, + BYD+ kS(B, + B;, + B,). : (10.12) 


Here D is the derivative symbol. Equation (10.11) gives 


(1 + k’S*)¢(z) = 1+ L7'Az — LF o(z) + Lo imck’ZS°N(2). (10.13) 
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Fic. 15e. Resonant length of a thin cylindrical antenna with flat ends and infinitesimal gap at 
different points. 
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Now the operator P inverse to 1 + k’S’ is easily found. It is 


(1+ kS°’)"f@ =Pf2eo=fe—k i} sin k(z — #)f(t) dt. (10.14) 


0 


This means that for an arbitrary function f(z) we have P(1 + k’S’)f(z) = 
(1 + k’S’)Pf(z) = f(z); this is easily verified. 
In particular, we have 


P.1 = cos kz, Pz = k™ sin kz. (10.15) 


Thus (10.13) gives 
o(z) = cos ke + L™’Ak™ sin kz — L™*PFo(z) + L7iwck?ZPS’N(z). (10.16) 
This is a transform of the basic integral equation (6.3). It is exact because in deriving it, 


we have not actually used the approximation (5.9). With (10.16) we are to associate 
the boundary conditions (6.4). The first of these is satisfied automatically. The others 


give 
cos kl, + L™Ak™ sin kl, — L7'T,,PFo(z) + L™ inck’ZT,,PS°’N(z) = 0, 
(10.17) 


cos kl, + L~'Ak™ sin kl, — L™'T,,PFo(z) + L™iack’?ZT,,PS’N(z) = 0. 


Here 7,, , T;, are substitution operators, meaning “put z = |, , or z = l, , finally.” 
We now eliminate A and Z from (10.16), (10.17); this gives for ¢(z) the equation 
| ¢(z) — coske + L™'PF¢(z) sin ke PS’N(2) 


| 
| 
| 





— cos kl, + L7'T;,PFo(z) sin kl, T,,PS’N(z) | = 0. (10.18) 





— cos kl, + L7'T,,PF¢(z) sin kl, T,,PS’N(z) 


Remembering that N(z) is known, the plan of solving by successive approximations 
is now obvious. We are to substitute @¢ = ¢, in the PF column, ¢, being some initial 
approximation, and solve, obtaining ¢, . Then ¢, is to be substituted in the PF column, 
and the equation solved, giving ¢; ; and so on. At any stage, we might get Z from (10.17) 
but it seems probable that a better value will be given by (6.23) or (6.25). 

Under the assumptions (7.1), (7.3) we have (5.9), and hence we can calculate PS*N(z). 


In general 


Pf(z) 


fa —k [ f(t) sin kf — tadt 


(0) cos kz + | Df(t)- cos k(z — t) dt (10.19) 


/o 


f(0) cos ke + k-*f’(0) sin kz + k™’ 4 D*f(t)-sin k(z — #) dt. 


Hence 


fad 


PS*N(@) =k | N(t)sin ke — 4) dt, (10.20) 
“0 
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Fra. 15d. Shape term in the reactance of a thin cylindrical half wave antenna terminated by 
spheroids of semi-axis b (center gap). 


and substituting from (5.9) we have 


lL<z2< —n:PS’N(z) 


—2k~*n~'[cos k(z + n) — cos kz], (10.21a) 


—n <2< 9 :PS’*N(z) = —2k-*n "(1 — cos kez), (10.21b) 


42S t 2 PS’N() = —2k*n~'[cos k(z — ») — cos kz]. (10.21¢) 

It is natural to take as first approximation ¢,(z), the function obtained by deleting 

the PF terms from (10.18). This ¢,(z) is precisely that given by (7.10). Then the second 
approximation ¢,(z) will be given by 


| ¢(z) — coske + L”'PF¢,(z) sin kz PS’N(z) | 


— cos kl, + L'T,,PF¢,(2) sin kl, T,,.PS’N(z) | = 0, (10.22) 


— cos kl, + L7'T,,PF¢,(2) sin kl, T,,PS’N(2) | 


and the higher approximations by similar formulae. 


APPENDIX: NOTES ON FIGURES 14 AND 15. 


Fig. 14. Relative current in a thin antenna. The following graphs are drawn from 
equations (7.10), which were obtained on the following assumptions: 
(i) the tangential electric field is constant over the gap; 
(iz) the antenna is very thin, but not necessarily cylindrical with flat ends; 
(777) when the gap is infinitesimal (as in Figs. 14 a, d, f) its length is still much 
greater than the infinitesimal radius of the antenna. 
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The ordinate ¢,(z) in each graph is the principal part of the relative current, J(z)/I(0). 
The abscissa z represents position on the antenna. The ends of the antenna are z = l, , 
z = I, , and the ends of the gap are z = —7n, z = 7, so that z = 0 is the center of the 
gap. Also 

k = 2x/yX, \ = wave length, 


21 = l, — l, = total length of antenna. 


In each figure graphs are drawn for various values of kl up to kl = z. 
The following points are of interest: 
1) When the gap is finite, the derivative of ¢,(z) is continuous; when the gap is 
infinitesimal, the derivative is discontinuous at the gap, unless kl = 32 or z. 
2) The curve for kl = 3} (that is, 21 = }4)) is always a single sine curve, whether 
the gap is finite or infinitesimal. 
3) Certain curves are not shown, because they go to infinity. This means that 


I(0) = O and the impedance is infinite’ (to this order of approximation). This 
occurs for kl = aw (that is, 21 = Xd) in Fig. 14a. This infinity disappears when 


the gap is widened in Fig. 14b. This does not mean that the widening of the gap 
eliminates infinite impedance; in Fig. 14b there is infinite impedance for kl = 
82/7, since this makes K(n) vanish in (7.11). However, when the gap is finite, it 
is questionable whether the impedance is correctly defined (for matching pur- 
poses) by Z = V/I(0). 

For infinitesimal eccentric gaps, we get infinite impedance for kl = 47/5 in 
Fig. 14d and kl = 27/3 in Fig. 14f. 

4) Fig. 14ce has a remarkable feature: the current is constant in the gap for kl = 7z. 

Fig. 15. Impedance of a thin antenna. These graphs are based on equations (8.10) 
and (9.1). 

Fig. 15a shows the resistance of a thin antenna plotted against kl. The gap is in- 
finitesimal and situated at the point of quadrisection. It is interesting to compare this 
curve with those given by King® and Schelkunoff (loc. cit.) for an antenna with central 
gap. The effect of moving the gap from the center to the point of quadrisection is to 
change the point of great or infinite resistance from kl = x to kl = 21/3. This comes 
from the factor sin kl, in the denominator in (8.10). Further, the resistance of a half- 
wave antenna (kl = $m) is changed by this shift of gap from about 70 to al out 140 ohms. 

Fig. 15b shows the reactance of a thin cylindrical antenna plotted against kl. As in 
the case of Fig. 15a, the gap is infinitesimal and at the point of quadrisection. Graphs 
are plotted for several values of the thickness parameter K = In (l/a), where 2l is the 
length of the antenna and a its radius. The reactance vanishes not only in the neighbour- 
hood of kl = 4 and kl = z, but also in the neighbourhood of kl = 27/3 = 2.09. - 

The resonant length of a thin cylindrical antenna (i.e. the length making the re- 
actance vanish) is a little less than half'a wave length. The shortening below the half 
wave length is a function of the position of the gap. This dependence is shown graphically 
in Fig. 15ce. 

Fig. 15d shows the effect of rounding the ends of a cylindrical antenna, as in Fig. 13. 
The antenna is half a wave length long and the gap is infinitesimal and at the center. 


Flat ends correspond to kb = 0 and a completely spheroidal antenna to kb = 3a. 





*L. V. King, Phil. Trans. Roy. Soc. (A) 236, 381-422 (1937). 
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THE RADIATION AND TRANSMISSION PROPERTIES OF A 
PAIR OF SEMI-INFINITE PARALLEL PLATES—I* 


BY 
ALBERT E. HEINS 
Carnegie Institute of Technology 


1. Introduction. We are concerned here with the following problem. A plane mono- 
chromatic electromagnetic wave is incident upon a pair of semi-infinite parallel metallic 
plates of zero thickness and perfect conductivity (see Fig. 1 for a side view). The edges 
of the plates are infinite straight lines which are parallel to the y axis of an zyz rec- 
tangular coordinate system. (The y axis is perpendicular to the plane of the paper in 
Fig. 1). The plates extend indefinitely in the direction of the positive z axis and are 
spaced a units apart. It is assumed, as in CHI’, that the electric field of the incident 
wave has only one component, namely the one which is parallel to the y axis. Since 
the incident electric field is independent of y, and the boundary conditions on the plates 
are fulfilled independently of y, no other components of the electric field will be excited. 
There will be two components of the magnetic field; these in turn may be derived from 
the single component of the electric field through the Maxwell equations. The angle 6, 
the direction of the propagation vector of the incident wave, is measured with respect 
to the positive z axis. 

We have just described the manner in which Fig. 1 is excited from free space. It 
is now necessary to indicate the mode of excitation which the parallel plate region can 
sustain. We assume that for z > 0,0 < x < a, the y component of the electric field is 
asymptotic to (p,e""* + p.e ‘“*) sin (xx/a). That is, the parallel plate region can sustain 
a mode which is consistent with the polarization which we consider here. It is to be 
understood that there are no other means of excitation in the finite part of the xyz 
space. x is the propagation constant in the parallel plate region and is equal to 








(k° — (w/a)*)'”*, where k = 2x/) and 2 is the free space wave length. In order that a 

(0,a) 
x 
A 
Ld 

(0,0) Zz 
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Fic. 1. 


*Received March 29, 1947. This paper was presented to the American Mathematical Society on 
Dec. 27, 1946. 

‘Carlson and Heins, The reflection of an electromagnetic plane wave by an infinite set of plates I, this 
Quarterly, 4, 313-329 (1947). Hereafter we shall refer to this as CHI. We employ here the same form of 
the Maxwell equations which were used in CHI. The time dependence is taken as exp (—ikct) where c 
is the velocity of light. 
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single mode propagate in the parallel plate region, it is necessary to assume further that 
1 < a/x < 1. The constant p, is the amplitude of the wave going to the right in the 
parallel plate region, while the p, is the amplitude of the wave going to the left in the 
same region. 

This structure may be viewed as a two dimensional antenna in the sense that it 
can receive and transmit energy. The problem can be formulated mathematically as a 
pair of simultaneous integral equations of the faltung type which are closely related to 
those of the Wiener-Hopf type.” The unknown functions in these integral equations 
are the surface current densities on the plates. The solution of these equations will give 
us the functional form of the current densities, their asymptotic form for z ~ 0°, z— @, 
as well as the relation between the amplitudes of the various waves. We shall divide 
this problem into two parts. In the first part we shall assume excitation from free space. 
Then the parallel plate region sustains the wave travelling to the right, since there are 
no obstacles in the parallel plate region which would give rise to a wave travelling to 
the left. In this case, we find the magnitude as well as the phase of the parallel plate 
wave which is travelling to the right. In the second part, we shall assume that the 
parallel plate region has been excited. Here we shall find the reflection coefficient, that 
is, the ratio of p,/p. . This second problem breaks down into a single integral equation 
due to the presence of a symmetry in the field components. The first problem we treat 
considers the parallel plate region as a receiving antenna, while the second one considers 
it as a transmitting antenna. We shall see that their properties are not completely 
independent. 

The formulation of the pair of simultaneous integral equations which we have just 
mentioned can be carried out by the same method employed by Carlson and Heins 
(CHI). An application of Green’s integral theorem in two dimensions with a free space 
Green’s function as a kernel gives us the y component of the electric field in terms of 
the surface current density on each plate./Thus if H,(z, z) is the y component of the 
electric field, and J,(z) and J,(z) the surface current densities on the lower and upper 


plates respectively, we have the following relation 


E(x, z) = E,” (2, 2) 


(1.1) 
| {T,(2’)H3’ [k(a? + (2 — 2’))'?] + Le’)Ho [k(x — a)? + (2 — 2’)*)'”7)} de’, 


2 
tS 
where H;” is the Hankel function of the first kind and E}"*(a, z) = exp [tk(x sin 6 + 
z cos §)]. The boundary conditions on EZ,(x, z) give us the simultaneous integral equa- 
tions. Indeed, since E,(z, z) is the component of the electric field which is tangent to 


the planes x = 0, z > 0 and x = a, z > 0, we have 


0 


Ey", 2) +7 | (Io(e)H@[k | 2 — 2’ |] 
: . (1.2a) 


+ 1,(2)HS? (ka + (2 — 2’)’)'”"]} de’, 


2Paley and Wiener, The Fourier transform in the complex domain, Am. Math. Society Colloquium 
Publication, 1934, ch. IV. Actually, the integral equations we are required to solve are singular cases of 
the Wiener-Hopf theory but they are still susceptible to Fourier techniques. 
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and 
zi" (a, z) + an {T,(2’)H5” [k(a? + (z om 2’)*)'”"] 
(1.2b) 
+ I,(2’)Ho’ [k | 2 — 2’ |]} de’ 
for z > 0. 
We can simplify these last equations by performing the arithmetical operations of 
addition and subtraction. Upon adding, we get 


0 = B"0,2) + Ba, +5 [Joe Hi tk |e - 2’ [ 
4 0 
(1.3a) 
+ HSK + @ — 29", 
while subtraction gives us immediately 
= R™(0, 2) — Ei™(a, 2) + ; [ Jie) {HO(k | 2 — 2 |] 
(1.3b) 


— Hy’ [k(a’ + @ — 2’)*Y'")}. 


Here J,(z) = Ip(z) + I,(z) and J,(z) = I,(z) — I,(z). In view of the z dependence of 
the kernels and the particular limits of the integrals in Eqs. (1.3a) and (1.3b), we have 
here two integral equations which may be solved rigorously with the Fourier transform 
in the complex domain. This implies, of course, that we seek solutions of appropriate 


growth, and the kernels possess the correct growth. We shall now show that such is 


indeed the case. 

2. The Fourier transform solution of equations (1.3a) and (1.3b). Let us now write 
Eqs. (1.3a) and (1.3b) in a form which makes them amenable to Fourier transform 
methods. We define £:"*(0, z) and E£;"‘(a, z), Jo(z) and J,(z) to be identically zero for 
z < 0. We further extend the Eqs. (1.3a) and (1.3b) for z < 0 to read 


ra 


éo(2) = § [Joe HO [k |e — 2 I + HOR + @ — 2") ae’, (2.18) 


i) 


; | Ji?) {He lk | 2 — 2 |] — HO tk? + ( — 2)?)'7)} de’, (2.1b) 


¢,(z) 


where ¢o(z) and ¢,(z) are defined to be identically zero for z > 0. Upon noting our as- 


sumptions on J,(z), J,(z), £,"*(0, z), E,"*(a, z), do(z) and ¢,(z) we have for all z 


(2) = Ei", 2) + Es(a, 2) + 5 | Joe’) HS? Ik | 2 — 2’ I) 
Jae (2.2a) 


+ Hs” [k(a® + (2 — 2’)’)'7)} de’, 
oi(z) = Ei(0, z) — Ej“(a, z) + 1 [ Ji(2’){ Ho” [k | z — 2’ |] 
Poe (2.2b) 
_ H$ [k(@?’ +(z- 2’))'”"}} dz’. 


” 
We assume as in CHI that k has a small positive imaginary part. 
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Some remarks on the growth of J,(z), J:(z), ¢o(z) and ¢,(z) as z becomes either 
positively or negatively infinite are now in order. With this information we can find 
the half planes of regularity of the Fourier transforms of the functions with which we 
have to work. It is to be noted, in light of the definitions we have imposed upon ¢,(z) 
and ¢,(z), that they are asymptotic to e ““*/z’” for z large and negative. This asymptotic 
form may be seen directly from the Hankel function. Thus the Fourier transforms of 
¢o(z) and @¢,(z) are 


nd 


. 0 , 
&,(w) = | e *”*b(z) dz and ®,(w) = / e *”"g,(z) dz 


and ®)(w) and ®,(w) are regular in an upper half plane ¥mw > —&mk. 
The transforms of the Hankel functions 1/4 H5”[k | z|] and i/4 Hj? [k(a® + 27])'? 
are well known. For example® 


Sf HS (Ka? + 2)7)e" de = EE — wy exp fi a | — wy?) 


and is regular in the strip —Smk < &Ymw < Smk. Furthermore the transforms of 
EE; (0, z) and E£;"*(a, z) are readily calculated since they have been annihilated for z < 0. 
Hence we have 

1 
i(w — k cos 8) 





[ ioe Ee z) dz =e 


“0 


and 


ak sin @ 


e 


[ e'* Ey” (a, 2) dz = iw — b cos @) 


Jo 
and these last two transforms are regular in the lower half plane §¥mw < Smk cos @. 
We observe, that thus far the transforms of ¢)(z),.¢,(z), the Hankel functions, £}"*(0, z) 
and F}"*(a, z) are regular in a common strip —Y¥mk < Y¥mw < Ymk cos 8@. 

We still have to discuss the growth properties of the surface current densities /,(z) 
and J,(z) for z > 0. Their dominant parts for z > 0 are terms of the type e'**. All other 
terms in the asymptotic forms of J,(z) and J,(z) approach zero more rapidly than these 
imaginary exponentials. Now a term of the type e**’ has the Fourier transform 


aco 


a 

Jo 
which is regular in some lower half plane bounded by a small but positive ordinate. 
It now follows that the Fourier transforms of ¢)(z), ¢,(z)Jo(z), Ji(z)E;"*(a, z), E{"* (0, z) 
and the Hankel functions are regular in the strip —SYmk < Ymw < Ymx« (or ¥mk cos 8). 
We are thus permitted to apply the Fourier transform to Eqs. (2.2a) and (2.2b) to get 


taksind . ta(k?—w?)t/9 
iid eo ote OS H,(w), (2.4a) 


(w—k cos 6) 2 w?)'/? 





(1 _ ” sanwond i (1 = gitt—erer . 
ake ror we wy H,(w), (2.4b) 








8The branch of (k? — w?)!/? is equal to k for w = 0. 
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where 
H,(w) = | e ‘"* J,(z) dz, H,(w) = [ e~*”* J,(z) dz. 
0 “0 


Equations (2.4a) and (2.4b) are now to be decomposed into two sets of equations, 
one of which will be analytic in the upper half plane {mw > —$mk while the other of 
them will be regular in the lower half plane ¥mw < mk or (¥mk cos 6). Let us first 
turn to Eq. (2.4a). The factor 


(k? — w*)"'?{1 + exp [ia(k? — w’)’”]} 
may be written as 
2 2) 1/2 -Dp2 _ , 21/2 @ 772 _ ,,2)1/2 _ KW) _ 
2(k w) exp E » (k w) | cos E (k w) ] "ra" K(w). 


Without indicating the precise form of K_(w) and K,(w) which we assume to be regular 
in the appropriate lower and upper half planes, we may proceed to the required de- 
composition of Eq. (2.4a). We have 


(1 + eK .(w) — K.(k cos 6)] 
i(w — k cos @) 





,(w)K.(w) — 
(2.5a) 





_ a+e"")K(k cos 0) , i 
7 i(w — k cos @) + 2 K_(w)Ho(w). 


The left side of Eq. (2.5a) is regular in the upper half plane ¥mw > —$mk, while the 
right side is regular in the lower half plane ¥mw < &mx (or ¥mk cos 6) and both sides 
are regular in a common strip. It follows then that each side is equal to an integral 
function €)(w), i.e. 

+ e***'"*)K .(w) — K.(k cos 6)] - 


b.(w)K.(w) — & pe = «(w), (2.6a) 





(1 + e***'"") K(k cos 8) 
i(w — k cos @) 





4 © K -(w)Ho(t) ali. (2.6b) 


In a similar fashion, we may decompose Eq. (2.4b). Let 


(ke — w) {1 — exp [ia(k’ = w’)'”*)} = e = Lv), 


where L_(w) and L,(w) are regular in the appropriate lower and upper half planes. 
We have upon repeating the argument for separation 
1 — eT L.(w) — L.(k cos 6)) _ 


b,(w)L,(w) — ( i—-tae = ¢,(w), (2.7a) 





° ae taksin@ % - 
5 L-(w)H,(w) + a ite — me . ®) = aw), (2.7b) 





where ¢,(w) is an integral function. 
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We are now compelled to indicate the precise forms of K_(w), K,(w), L_(w) and 
L.,(w) if we are to evaluate the integral functions €)(w) and ¢,(w) and thereby find 
H,(w) and H,(w). In much the same manner which was indicated in CHI one finds that 


2 (w-—k de a lb 1/2 
K_(w) = = - i exp E (k? — w’)'” arc tan (rte) 


Ps = =)” taw | —saw/ x (2n+1) 
+ xo(e) | IT (1 San+t) +zm+ph 


n=l 
is regular in the appropriate lower half plane.* x,(w) is an integral function which has 
been introduced into the product decomposition of K(w) and is to be chosen such that 
K_(w) is of algebraic growth for | w|-—© and $mw in the correct lower half plane. 
Similarly 


1 _ Qw+x) E 2 tye, fi=sy" 
K,(w) = (k +p w)'”? exp 7 (k w) are tan a ee 





= (| Il (1 — — (chy) =m ___taw | on 
aa 3 (2n + 1)° r(2n + 1) ; 


n=1 
In order to determine x.(w), we simply calculate the asymptotic form of K_(w) 
and K,(w) as |w|— ©, Ymw in the appropriate half plane and choose x,(w) so that 
K_(w) and K,(w) will be of algebraic growth.° Let us first study K_(w). If we observe 
that the parameter ak may be neglected in the infinite product as |w|— ©, ¥mw < 
Ymk, we have that K_(w) is asymptotic to 


1/2 oe | tw 1 = )| T[ daw | -seoretmen 
w exp | i log ( ke + xo(w) I] oT ae + bf 


= ew? Lliaw/2m) . | = ( ee =) | 
= cw T(iaw/n) exp | — 1 9 + 9 log ke + xo(w) 





(2.8) 


where c is a constant whose precise form does not interest us and y is the Euler-Masche- 
roni constant. We may now apply the Stirling formula to (2.8) to obtain that K_(w) 
and K,(w) are of algebraic growth for | w| — ©, $mw in the appropriate half planes, if 





aw 


xo(w) = > | -3 "y= log = - ir]. 
With x,(w) so chosen, K_(w) is asymptotic to w'” for |w|— © and §mw in the 
appropriate lower half plane, while K,(w) is asymptotic to w'” for | w| — © and §mw 
in the upper half plane ¥mw > —Smk. 

We are now in a position to determine the integral function «,(w). Let us note in 
Eq. (2.6b) that H,(w) is the unilateral Fourier transform of J,(z) + I,(z) = Jo(z). As 
such, since J,(z) has appropriate growth for z large and positive, and since it is integrable 
over any finite interval of z including the origin, H,(w) possesses the property that 


‘Henceforth principal determinations of inverse trigonometric functions and logarithms are under- 


stood. 
5J. S. Schwinger, Theory of guided waves, Radiation Laboratory publication, forthcoming. 
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it approaches zero for |w|— ©, ¥mw < Ym« or Ymk cos 6. If we let | w| > &, w 
in the correct lower half plane we see from Eq. (2.6b) that ¢(w) = O(w~*~’””), a > 0. 
On the other hand ®,(w) approaches zero for |w|— © SYmw > —S&¥mk because #,(w) 
is the unilateral Fourier transform of a function defined for negative z, is integrable 
over any finite negative range of z including the origin and possesses appropriate growth 
for z large and negative. Hence for |w|— ©, §mw > —S&mk, we see from Eq. (2.6a) 
that ¢,(w) = O(w'?~*), B > 0. It follows, therefore, by a theorem of Liouville, that 
€o(w) is a polynomial of degree less than minus one half, and hence identically zero. 


We have finally 


2(1 + e***'"")K .(k cos 6) 
(w — k cos 0=)K_(w) ’ (2.9) 





H,(w) = 


the Fourier transform of Jo(z). 
We can obtain some information regarding J,(z) for z — 0° from H,(w) as | w| —- 
o, ¥mw < Ym« or ¥mk cos 6. For now we have that 


H,(w) = O(w~’”’) 


and this tells us immediately that 
Jo(w) = Of") 
for z — 0°. This verifies the integrability of J,(z) for finite and positive z. 


We now turn to the determination of H,(w) and for this we consider Eq. (2.4b). 
We note that 


a = — i(k? — w)”? exp E (k? ae uty | sin [g (k? = wt) | 





where now 





° 1/2 
L.(w) = —ia exp E (k? — w*)'” arc tan (i + “) 


+00) fill - )" +g 


and 





1 ee za \—— ( _ w\" 
a ae exp E (k w’)’” are tan ye 


-w06] f[(- @))” 


Here we have the L_(w) which is regular in the lower half plane ¥mw < $mk and the 
L,(w) which is regular in the upper half plane ¥mw > —Smk. Again we choose x,(w), 
an integral function, such that L_(w) and L,(w) will behave algebraically as | w|— 
and $mw in either of the appropriate half planes just described. We proceed as we 
did above to determine x,(w). For |w| —~ ©, S¥mw < &mk, L_(w) is asymptotic to 
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exp | xu) + = log — a | Il a+ taw —itaw/2rn 


=< *s. 


_ 2x exp [x:(w) — tawy/2e + {iaw/2r} log (—2w/k)] 
pa iawl (taw/2r) 





Upon applying Stirling’s expansion theorem we obtain immediately 


A exp [x:(w) + {iaw/2r}{1 — y + log (— 2w/k)})( Sa) sil 


iaw 

Hence, if we choose 
7aw ix 
= — l 

L_(w) will have algebraic growth in the lower half plane §¥mw < &Ymk for |w|— 
and will be asymptotic to w~’””. If we repeat the argument in the upper half plane, 
Smw > — Ymk for the term L,(w) we find that the same x,(w) will render L,(w) 
algebraic in growth for | w|— © and now L,(w) will be asymptotic to w'”’ for w > o. 
Finally, reasoning as we did for Eqs. (2.6a) and (2.6b) we find that ¢,(w) is O(w*'), 
a, < —}for|w|— ©, ¥mw < $m and is O(w"), 8, < } for | w| — ©, ¥mw > 
—¥mk. Applying Liouville’s theorem once again, we find that ¢,(w) is identically zero. 
Thus we now have 


(1 — e****'"*)L..(k cos 6) 
(w — k cos 6)L_(w) 





H,(w) = 


the Fourier transform of J,(z) — J,(z) = J, (2). 

From H,(w) and H,(w) we can obtain the Fourier transforms of J,(z) and J,(z) by 
a simple addition and subtraction. We can also see how J,(z) behaves for z — 0°. Since 
H,(w) is now O(w*””) for |w| > ©, ¥mw < Ymk, J,(z) = O(2"'”’) for z — 0* and 
this verifies the integrability of J,(z) for finite and positive z. In closing, we note that 
the precise forms of J,(z) and J,(z) are of no interest to us since we are only interested 
in the far fields. 

3. The calculation of the far fields. In order to calculate the far fields we first 
express Eq. (1.3) by a Fourier integral representation. We have 


E(x, z) = exp [ik(x sin 6 + z cos @)] 
7 tw2z/7,2 2\-1/2 i[z| (k29—w2)2/* 
+z f ee — wy [{Holw) + Hi(w)}e 
8r c 


+ {Ho(w) — Hi(w)je' #9") dw 


where C is a path of integration drawn within the strip of regularity of all the Fourier 
transforms which appear in the above integral. The path is closed either above or below 
depending upon whether z > 0 or z < 0. Care must be taken in closing the path so that 
it does not intersect the branch cuts which are introduced due to the presence of the 
branch points k and —k. The dominant terms arise from the residues due to the two 
poles k cos @ and x. Furthermore, since contributions from other poles or branch points 
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give rise to terms which are small compared to the terms arising from the poles k cos @ 
and « for |z| —> ©, we need only calculate these dominant effects, at least insofar as 
we are concerned with the far field.° 

There are four separate regions of interest (a) z < 0, —~ <4 < @ (b)z>0, 
x >a(c)z> 0,2 < Oand (d) z > 0,0 < xz < a. Let us consider region (a). For z < 
0, E,(x, y) is asymptotic to exp [ik(z cos @ + z sin @)], as it should. For region (b), 
z> 0, z > a, E,(z, z) has no term comparable in magnitude with the plane wave term. 
Thus region (b) is the region of the geometrical shadow. For region.(c) z > 0, x < 0, 
E,(x, 2) is asymptotic to 


2i exp [zkz cos 6] sin [kz sin 6]. 


Thus for z < 0, z > 0, the lower plate acts as a perfect reflector. 

Region (d) is the interesting one. We now have a means of finding the amplitude of 
the transmitted wave guide mode. We know that in this region E,(z, z) is asymptotic 
to e'“* sin rx/a for z > 0,0 < x < a. On the other hand when we evaluate the integral 
and take out its dominant terms, we find that we are left with 


K(k cos 6) im we eft in aii 
x — k cos 0... K_(w) erie 





ae au (1 + ces 
T 


The amplitude of the transmitted wave is the coefficient of the factor e** sin rx/a. It 
may be simplified if we now take k to be real. In the first place 





lim 5 = 2[x°(k — «)a~*]'”? exp | io — xo(x) — are tan (; : ‘)"] 
wv K_(w) 2 1 ) = 


where 





— P Ka Ka 
0, = —- are sin - = | 
>| a[(2n + 1)? — 1}? 3=(2n + -1) 
while 
K.(k cos @) 


_ a{k(1 + cos 6)(k? cos? @ — «’)]'” exp [70. + xo(k cos 6) — (iaké sin 6)/2r] 
2r(k cos 06 + «)[cos {(ak sin 6)/2}]'” 





and 





bi > | are si ak cos 6 _ _ak cos 0 | 
[Qn + 1)? = (ak sin 6°)? 72n + IT 


n=1 


Thus, the transmission coefficient is 





kr(1 + cos ) cos (ak/2 sin af tV+a(kceos0—x)/4 
1 a’(k* cos’ @ — « )(k + x) . » OS OSr 


6—. T. Copson, Oxford Quart. Math. 17, 19-34 (1946). There is a detailed discussion in this paper on 
the choice of a path similar to C. 
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where 


- k= 1) ak (, _ 8), 
¥ = 0, + 0, + arctan (L—) “—" 1 — ~)sin 6 





+ = = [y — 3 — log (x/ak)]. 





The square of the absolute magnitude of the amplitude of the transmission coefficient 
is proportional to the power gain, as a direct consequence of the Lorentz reciprocity 
theorem. Thus insofar as the angular variation is concerned, the radiation pattern is 


(1 + cos 6) cos (ak/2 sin @)e*"”’” 
(k? cos’ @ — x’) ; 





We have thus found how the parallel plates act as a receiving antenna. It is to be noted 
that the radiation pattern arises from the excitation of the parallel plates for z > 0, 
0 < z < a. The reciprocity theorem has enabled us to give a partial solution of the 
second part of the problem. The reflection coefficient which we have described earlier 


has yet to be calculated. 
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A PRACTICAL METHOD FOR SOLVING HILL’S EQUATION* 


BY 
L. BRILLOUIN 
Harvard University 


1. Introduction. The differential equation known as Mathieu-Hill’s equation can 
be written 


y”’ + J(z)y = 0, (1) 


where J is a periodic function of x. The period is usually taken equal to x for historical 
reasons. The first equation of this type was discovered by Mathieu in connection with 
the problem of vibrations within an elliptic boundary, when 


J(x) = n+ y cos 2x 
has a period z. 
Floquet proved that the general solution of Eq. (1) can be written 
y = Dye" &(xz) + Die“ 8(—2z), (2) 


where ® is a periodic function, with the same period w as J(x). This general solution 
contains two terms with the exponents +yu. Floquet’s theorem can be expressed in a 
slightly different way. Let us consider the term with +n: 


f(z) = e*®(z). (3) 
The condition on f(z) is 
f(z + nm) = e"f(x) = & f(x) 

with & = ¢”". 
The general solution is 

y = D,f(x) + D.f(—2z). (4) 
We shall look for a solution in an interval of length 7 and use condition (3) to extend 
the solution from — © to +. This means that we shall have to meet some boundary 
conditions in order to match the solutions in two consecutive intervals. These matching 
conditions will be essential in fixing the value of u. 

The method developed in this paper is based upon these general considerations and 
shall be explained more completely in Sec. 2. It differs completely from the classical 
method, as found in most textbooks.’ The standard procedure is to expand J in Fourier 
series 


J _ =. a (5) 
and to look for a similar expansion for the unknown periodic function ® 


@= > be™. (6) 


n 








*Received Feb. 4, 1948. 
1Whittaker and Watson, Modern analysis, Cambridge University Press, 4th edition, 1927, p. 414. 
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Equation (1) results in an infinite system of simultaneous linear homogeneous equations 
for the unknown 6, s. A non-trivial solution can only be obtained if the corresponding 
infinite determinant is zero. This last condition is used to determine the value of the 
exponent u. 

Whittaker was able to discuss this condition in the case when the series of the co- 
efficients @, in expansion (5) is absolutely convergent, and obtained a formula 


sin? (z in) = —gsinh? (z u) = A,(0) sin’ (z ) (7) 


where A, (0) is another infinite determinant with the following coefficients 
A, (0) — | Bus , } = E, 
6... (8) 


~ 0 — 4m” allt ae 








| 


This result is not very encouraging. First, the condition of absolute convergence for 
the series of the coefficients of 6, is a very restrictive one. Second, we still have to com- 
pute an infinite determinant, and the computation proves very difficult unless the 6, 


terms decrease very rapidly when n increases. 
This infinite determinant takes on infinite value whenever 


Oo = 2n ninteger (9) 


These 6, values correspond to double poles of the determinant, since both rows m = +n 
obtain infinite terms. These double poles are canceled out in formula (7) by the double 
zeros of sin’ (76,’"/2) and do not correspond to any singular values of u. 

Altogether, the method of Fourier expansion is not very practical, and leads to 


complicated computations. 
The method presented in this paper does not involve the restrictions of the classical 


method, and leads to a practical solution of Hill’s equation, even in such exceptional 
cases as periodic functions J(x) containing discontinuities or 6 functions. 
2. Principle of the method. Let us consider a differential equation 


y”’ + F(z)y = 0 (10) 


with a given function F(x). We may find two independent solutions u and v and obtain 
the general solution 
y = Au+ Bo (11) 
containing two constants A and B. We obviously have 
uo’ = vu” = Pu, 


hence 
uv’ — vu’ = C, 


and a suitable normalization of u and v is used to make the constant C unity: 
(12) 


u’ — vu’ = 1 
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THE PERIODIC J(x) FUNCTION 
Fia. 1 
We now want to discuss Hill’s equation 
y’ + J(xz)y = 0 (13) 


with a periodic function J of period defined in the following way (see Fig. 1): 
F(x) —7nr/2<2< 2/2, 
F(x — nx) n/2< 2x < 3n/2, (14) 
F(x — nr) nx — 2/2 <2 < ne + 2/2. 


Floquet’s theorem assures us of the existence of two independent solutions y, , y. char- 
acterized by the following properties (Eq. 5): 


yi(x + er) =e” y,(x) = Ey (2), 


&~'y,(z), : (15) 


y(x + r) = e“"y,(z) 


ur 


g=e 


Let us consider y,; and discuss a practical method for obtaining u (or &). In the first 
interval (—x/2 < x < 2/2) the function y, may be represented by a formula (11), 
with a set of A and B coefficients. In the second interval (x/2 < x < 32/2) the co- 
efficients will be £A and &B, according to (15). We must now write the continuity con- 
ditions for y, and y{ across the border 2/2: 


Au, + Bo, = Au, + EB, , 
Au} + Bui = tAus + EBoi , (16) 


u, = u(r/2), Uz = u(—7/2), v, = v(r/2), Vo = v(—7/2). 
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We have obtained a set of simultaneous equations for the two unknowns A and B. 


This can be solved only when the determinant is zero: 


uU— EUs = fv, 
D= 
ui — kus vi — £v3 (17) 


= 2 + tulv, + uy, — uvi — uri) +1 = 0, 


where we have used Eq. (12) at both points +2/2. Equation (17) in é fixes the Floquet 
coefficients £ and &*. The product of the two roots is unity and their sum is given by 


2 cosh ur = §€+ €° = —ulve — ubv, + uns + ugvt; (18) 
hence 


4 sinh’ ( *) = 2cosh yr — 2 = —(u, — w)(vi — vi) + (ui — wi), — v2) = (19) 


with the help of Eq. (12). Once the wu exponent is obtained, the coefficients A and B 
result from (16) and the solution of Hill’s equation (13) is achieved. 
A very important case is obtained when 


F(x), J(z) = even. (20) 
One may choose correspondingly for u, v an even and an odd function: P 
, , 
u(r) = u(—2), Ue = U,, u= —uU, 
(21) 
v(x) = —v(—2), ve = —v,, ve = vf 
Eq. (18) then reads 
cosh ur = Uyvl + usr, (22) 
from which we obtain 
sinh’ (um/2) = ufr, , 23) 


a formula that will be found useful for a comparison with Whittaker’s theory of Hill’s 
equation. 

The whole method is a generalization of the discussion that was previously given 
for the special example of a rectangular J(x) function.* One advantage of the method 
is that it works with periodic J(x) functions exhibiting a finite number of discontinuities. 

3. Some special examples. The case of a rectangular J(x) function can be easily 
investigated along these lines, and the results agree completely with those of a previous 
discussion’ where a slightly different method was followed. 

Let us consider the case of a parabolic function 


F=a-— b’x’. (24) 
*L. Brillouin, Wave propagation in periodic structures, McGraw-Hill, New York, 1946, Ch. VIII, 


pp. 180-186, and Ch. IX, pp. 218-226. 
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We may find the solutions as power series expansions 


“u= 1 + uz’ + coer Uo” -°*, 


(25) 
vp = aot vet tees dong tt --- 
Substituting in Eq. (10) we obtain the recurrence formulas 
2u, +a = 0, (2n + 2)(2n + 1)ttans2 + Atle, — D'rnn-2 = O, 
(26) 


6v, +a = 0, (2n + 3)(2n + 2)von+g + Boner — D'Von-1 = 0, 


and the u, v functions satisfy the normalizing condition (12). The following special 


cases may be of interest 


a= —b, u=e” 
‘ (27) 
vy = eI (x) with I(x) = / e~”** dz, 
0 
a= —3b, u =e”? + Qbr eI (z), 
(28) 
y = xe’™” 2 


as may be checked by direct computation. 

In all these cases, there is no discontinuity of the function F on the limits +2/2 of 
the interval, and the curve on Fig. 1 is a continuous curve with a discontinuous de- 
rivative at +7/2. 

The corresponding Hill’s problem is immediately solved with the help of Equation 
(22) or (23). For instance, the cases indicated above under (27) and (28) yield: 


a= —b, sinh? (. *) = uUtv; = b : evry?) (29) 
( T r b T r b wT 

> le ee 2 2 q = a r*/a7i * Ee 

a 3b, sinh C *) b 4 + bre (?) +b 4° (%) (30) 


These results could not have been obtained by any other method of solution. 
Another example can be solved with the help of Bessel functions, which satisfy the 


equation 


a ( d ) eo 
. dz \ dz y) +@ n)y = 0, (31) 
y = AJ,(z) + BJ_,(2), n non integer (32) 
Taking a new variable 
x = logz, z=e', 


we obtain an equation of type (10): 


2 


d 
fy + Ct = ny = 0. 3) 
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Thus, 
y = AJ,(e*) + BJ_,(e*) (34) 
and we have a solution corresponding to an unsymmetrical 
F = e* — n° (35) 


for which our general formulas (18), (19) should be used in connection with the corre- 
sponding Hill’s equation. 

4. Solution with the B. W. K. method. The point of departure of our method is a 
solution of an equation of type (10). An approximate solution can be found with the 
B. W. K. procedure,’ if F(x) is exhibiting only small variations about a large average 
value. A more pretvise statement of the conditions involved will result from the fol- 


lowing discussion. We rewrite Eq. (10) as follows: 


y’ + G(x)y = 0, F(x) = G*(z). . (36) 
We now consider a function 
y=GQ'’e', S = [ G dx, (37) 
which yields 
y’ _ 3 () _1G@"_@ (38) 
y 4\G 2G 


The function y represents an approximate solution of Eq. (36) if the first two terms in 
(38) are negligible in comparison to the last one. This is the case if 


Ys Gg” 


a ™ €, we, e <1, (39) 
7 A, 


and terms in e’ are neglected, while ¢ terms are retained. The second condition (39) 1s 
very restrictive, however, since it allows only for variations of G’ of the order of €. 
The function (37) becomes a rigorous solution of Eq. (36) when 


‘ v/\2 ves 

3 ( ) 7 1 G aie (40) 

4\G 2G 

A 
y= — —e (41) 

: (a + x) 
Solutions u, v normalized in accordance with (12) are easily found: 

u = (2G) '’e"**, v = (2iG)*e*** (42) 


An interesting example is shown in Fig. 2; it corresponds to the even function 


G=-: f 2 (43) 
(a+ |x }) 


*The initials B. W. K. refer to the three authors L. Brillouin, Journ. de Phys. 7, 353 (1926); G. 
Wentzel, Zts. f. Phys. 38, 518 (1926); H. A. Kramers, Zts. f. Phys. 39, 828 (1926). 

For a complete discussion and more references, see E. C. Kemble, The fundamental principles of 
quantum mechanics, McGraw-Hill, 1937, Ch. III, p. 91-112. 
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with a discontinuity of the derivative at z = 0. Here, 





; sain 1 1 Ax 
— S(—7) = S ; = A = —-/ —_ — — = ——_-, 
S(—z) (x) [ G dz i(- = 1) da+s 


“0 


For the odd solution v we simply take 


v=G'’ sin S = en ! sin S, 
(44) 
v’ = -5 G ; *“G’ sin S ae G 2 cos S. 
Both v and v’ are continuous at xz = 0 since G’(+0) = —G’(—0). The even solution u 
can be written as follows: 
u=G *l cos S + K sin S], 
(45) 
; l y~3a/5 Y r ° ¥ ° Y r 
= —- GG [cos S + K sin S] + G'”’[—sin S + K cos S], 


y= 
2 


where the + sign must be taken for 0 < x < 2/2 and the minus sign for —/2 < x < 0. 

The function u is continuous at the origin, and u’ becomes continuous when it is made 

to vanish at the origin. Hence, 

cK -&@_ _@ 
2G°(0) A’ 

Finally, after some reductions, one finds 


= % 12 | | cos S + 4 sin s|, 


d 


(46) 





;-1|_l*el 2g A ++eet+/s).. | 
[tek oo . A(a + | x |) ain 3} 


It is easily verified that the coefficients have been chosen correctly so as to satisfy the 


normalizing conditions (12). 
With these solutions u, v we may now compute the Floquet coefficient wu, with the 


help of Eq. (23). We find 
A’? +a(a+|z}) 


= eee Se " ee 
sinh (u *) = uv; = - r sin S, cos S, — sin’ 8, re : 





where we must take 


An/2 
z= S, = —“7/-_ 
whe a(a + 2/2) 
sinh’* (, *) = “7 sin S, cos S, — Ate pane sin’ S, . (47) 
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A simple check can be made on this formula. Keeping the same A value, we may replace 
a by 
es Tv 
7 =~ 2 
thus interchanging the position of the angular maxima and minima of the curve (see 
Fig. 2). Our formula is not affected by this change. 
5. Successive approximations. Starting from the solutions obtained in the preceding 
sections, we may develop a method of successive approximations. Let us consider Eq. 
(10) with a function F and assume that it can be approximately represented by a G’ 


function of the type (41). 
A 2 
= —————,, F =G@ —- H(z), small. 

G (a + a) G eH (x) esma (48) 
If a single function G does not yield a good enough approximation over the whole 
interval —7/2, +72/2, it may be convenient to divide this interval into two or more 
partial intervals using different G functions, and to join solutions at the boundaries, as 
shown on Figs. 2 or 3. 

Next we use an expansion 


Y=Yrryteyte. 
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Grouping terms with the same power of e, we obtain 


ys’ + Gyo = 0, (49a) 
yt + @y, = Hy, (49b) 
yz’ + Gy. = Hy, . (49¢) 


The solutions of (49a) are the wu, and v, obtained in (42). Thus, the solution of (49b) reads 


z z z 
Y, = Uo [ Y,vh dx — vo [ Y,us, dx with Y, = [ Hyp dx (50) 
“0 bt) “0 
and the further approximations can be obtained in a similar way. 
As an example, let us assume an even function F which can be approximated with 
the even G function (43). The zero order approximation is represented by up and v of 
equations (44), (45). We thus have, — 


Ur = | Hud, Vo= {Hoy de. (51) 
Jo 0 
H being even, U, is odd and V, is even. Next, we obtain 
U, = U [ Uyvi, dx — v% [ U,us dx (even), 
“0 “0 
(52) 
Vv; = Uo [ Vv, dx — 0 [ Vous dx (odd), 
“0 “0 
and the first order solutions are 
U= Uy + eu t:::: (even), 
° (53) 
v = Uo oe €v; cae os (odd). 


These functions being automatically normalized according to (12), as a direct check 
easily shows. 

A solution of Eq. (10) can thus be obtained step by step to any desired degree of 
approximation, and may be used to the solution of Hill’s equation as shown in Sec. 2. 

The same method could be applied to any other known solution of Eq. (10). We 
used the G* functions in the zero order approximation, but any other known solution 
would do just as well. 

6. Hill’s equation containing delta function. Let us consider our fundamental Eq. (10) 
and assume the function F to be a delta function: 


(0 «#0 


F = Baz), 8 =4 , | o@ar=1 (54) 
| @ x=0 Ais 
Equation (10) can be readily integrated: 
(0, (x # 0), 
y 34 . 


‘|Bawye) (=O), 
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Thus, 


(a, (x < 0) 
= (55) 


la+By(0), (x>0) 


and we obtain our even u or odd v solutions: 


u=1+ 5B | x | (even), 
(56) 


v=2 (odd). 


These u, v functions are normalized according to (12). The solution of the corresponding 
Hill’s equation is obtained from Eq. (23): 


se T eu ] T 
sinh (1 "| = (uv) es/2 = (3 Bz) ae = 42 (57) 
This result can be checked by different methods. The case of a rectangular J function 
was discussed in the author’s book,’ p. 181, assuming 
J = (, + l, = m) (58) 
—x> (0 <2 < l,) 


The solution was given by equation (44.12), loc. cit. p. 181: 


1 , = , a 
cosh wr = cosh x,/, cosh x2l, + 3 (x “+ xs) sinh x,/,; sinh xl. . (59) 
“= \X Xi 


2 


We now take 


x, = 0, ~%—7e, 1-0, lx; = B, 
and obtain 
3 
cosh wr = 1 + : XehixXel, = 1 + B i, (60) 
2 2 
where 1, = w when /, = 0. This checks with our previous result (57). Other types of 


J functions can be used and lead to similar results. Such a problem is completely outside 
the reach of the Fourier series expansion method. 

7. Comparison between the present method and the classical one. We discussed 
in See. 1 the classical method of solution, and underlined its limitations. The Fourier 
expansion of the periodic function J must be such that the series of the Fourier co- 
efficients be absolutely convergent. This rules out functions with discontinuities, whose 
Fourier coefficients decrease as slowly as 1/n, but the classical method should apply 
to a continuous function J with discontinuous derivative. This is the case for the 
problem discussed in Sec. 4, Eq. (43): 


fe ee (61) 
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Within a period x this function oscillates between A?/a* and A?/((a + 2/2)*). The 
corresponding periodic function J can be analysed in Fourier series and Whittaker’s 
solution obtained. The question now is to compare solutions computed one way or 


the other. 
The comparison is easier when the variation of the function is small; hence we assume _ 


A = Ba’, a >: (62) 


and compute expansions with respect to the small quantity 


T 
€=5, <1: 


Let us start with the solution obtained in Sec. 4. It contains the quantity S, (Eq. 47). 


Aw/fB . BefB ef. ¢ zy... 
s, = (Ate = Bel = BE (1 z+ (z ) (63) 





we shall keep terms up to 2/2a only, since the computation of second order terms proves 
rather cumbersome. The solution is given in Eq. (47), and contains 


sin S, = sin (xB/2) — B (x*/4a) cos (xB/2) + ++: , 
cos S, = cos (rB/2) + B (x*/4a) sin (wB/2) + ---. 


Using these expansions in Eq. (47), we note a first term in (x/2A) sin S, cos S, that can 
be dropped, since A is of the order of a’. We thus are left with 


a(a + #/2)) 


° TT - 2 
sinh? (. 4 = —sin’ S, (1 + Ba! 


2 


and the term in parenthesis is again of second order. Thus, 


2 
sinh’ (1 =) = —sin’ S, = —sin’ (rB/2) + 2B 7 sin (xB/2) cos (rB/2) 
(64) 


= ain? R* wr 
sin B>3 +B 4a sin Br 


Terms in 1/a’ could be computed here without much trouble, but they lead to serious 
complications with the Whittaker’s formula, that we want to discuss now. 

We first compute the Fourier coefficients of the periodic function J derived from F 
according to Eq. (14): 


J= > oe" (65) 
Ne obtain 
2 «/2 —it2nz 2 a/2 —t2nz 
0, =.0..= = | —* __ dz = ge 2 f an Rowen 
w J-2/2 (1 + | x |/a) 7 /o (1 + 2/a) 
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The result is 


_ 2B" -( z)*| = ( _f 10 (x). .) 
9 aap +? = B\! a" 3 a 








(66) 
= _ 2B ae . on 
oi One ran” 1 (—1)" + (—1) s | 
Comparing (66) and (63), we note that to the first order 
S, = = oy” (67) 


9 8 

We also see that @, decreases as 1/n’ as needed, and that it is of the first order in 1/a 
Whittaker’s formula (7) reads 

sinh? (. r) = —A,(0) sin’ (z ai), (68) 

We shall prove easily that A,(0) is practically unity and equation (67) shows that Whit- 


taker’s formula (68) checks completely with our solution (64). 
Let us now discuss the infinite determinant (8): 


4,0) = |8.,|, B8..~= 1, 
die (69) 


%— 4m?’ (m ¥ p). 


Buy = 

Diagonal elements are all equal to unity while non-diagonal terms are proportional to 

6,,-» ; hence, according to (66), these non-diagonal terms are all very small, of the order 1/a. 

Such a determinant can be computed in the following way: we first take the product 

of the diagonal elements, that is 1. Next we take all the diagonal elements but two, 

namely n, n and m, m, which we replace by the non-diagonal terms B,,, and B,,, , then 

we take all the diagonals but three (n, n; m, m; p, p) which we replace by Bin, Buy ; 
Bim , ete. Thus, we obtain 


A — | Bus | _ 1 a a os + | a a a ml en 
(70) 


nFmMA PFN eecevce 


The rule is obvious, and the expansion is ordered with respect to powers of 1/a, with 
no term in 1/a and terms in 1/a’ , 1/a*® --- . We decided not to use terms in 1/a’ our 
expansions, hence our determinant is practically unity. Some difficulty may occur 
when 6, — 4m” becomes very small (of the order of 1/a), when the determinant becomes 
very large. We already noticed in Sec. 1 the inadequacy of Whittaker’s formula (68) 
near the poles of the determinant. There is a compensation, when 6, = 4m’, between 
the determinant having a double pole and the sin’ (4/2 (@)'””) a double zero. Our formula 
(64) does not exhibit any such trouble. 

Otherwise, both methods check completely. It is hoped that the general method 
developed in this paper will be found useful for practical discussion of many problems 


reducing to Hill’s equation. 
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—NOTES— 


THE SOLUTION OF NATURAL FREQUENCY EQUATIONS 
BY RELAXATION METHODS* 


By J. L. B. COOPER (Birkbeck College, London) 


The application of Dr. R. V. Southwell’s relaxation method’’’ to the solution of 
natural frequency equations for systems with a finite number of degrees of freedom is 
discussed in this article. The first object is to discuss the conditions under which the 
procedure converges. It is shown that two criteria can be given, such that if the one is 
satisfied, the procedure converges to the highest mode, and, if the other is satisfied, to 
the lowest. This is of more than purely theoretical interest: for it is shown that the 
possibility of finding the highest mode waiiied can be used to simplify considerably the 
process of finding further modes. 

1. The criteria of convergence. The equations to be solved will be written in the form 


D ant, = yin. (ry = 1,2---n) (1) 


where a,, and b,, are constants, | a,,| and | b,, | are symmetrical and the values of p* 
for which a solution is possible must be determined. We write 


= A(x) = > a,,2,2. , B = B(x) = >> d,.2,2, , 





0A 0B A 
a, = 3.) = % et 








Both A and B are positive definite forms. In general, \ denotes a variable function 
of the variables x, , but the letter \ with subscripts will be used to denote constant values 
corresponding to specified selections of the variables z. 

Suppose that x, is changed to x, + Az, . The change in d is then 


_A+4A_ A_ BAA — AAB 
‘= 24+ AaB B~ BB+ AB) (2) 


Now, 
BAA — AAB = B{A,Ax, 4 a,,(Ax,)°} - A{B,Ax, + by, (Aa,)*} (3) 


In the relaxation method, a first guess is made to select an initial set of values z, 
and the corresponding value of \, say ), , is calculated. Then, in the basic application 
of the method, one of the x, is adjusted to make A, — ,B, zero: the x, corresponding 
to the largest A, — ,B, is chosen. (In theory this is the best possible single adjustment 
that may be made; in practice it is rarely used as more useful operations can always be 


performed, according to the skill of the computer, who should always be thinking as 


*Received Aug. 29, 1947. 
1Southwell, Relaxation Methods, Oxford University Press, 1940. 
2Pellew and Southwell, Proc. Roy. Soc. (A) 175, 262-290 (1940). 
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far ahead as possible and therefore using his judgement to modify each adjustment he 
performs so that it will ultimately have the best possible effect when taken in con- 
junction with the adjustments which follow). The process is repeated a few times; then 
the value of \, corresponding to the values x so obtained, is calculated and used as the 
new value of \, in the next stage of the process. If the initial values x are correctly chosen, 
the values of \, found by these successive steps converge to the lowest value of p’ for 
which the equations (1) have the solution. We shall now proceed to give a rigorous 
proof of this last statement, and also to show how the process can, just as easily, be 
made to converge to the highest mode. 

The value of Az, , the change in x, , is determined by the condition that A, — A,B, 
vanishes when x, + Az, is substituted for x, . Hence 


A, — \.B, + 2 (Ger = Aibix) Ary = 0. (4) 
On substituting this value of Az, in (3), we get 


2(BA, — AB,(ay — xb 
(a2) Bau _ Abx) = ( k , BNC 1 | 





BAA — AAB 


ll 


(5) 





— Blan (0 sata Ai dix) -. (A — mo) ba ae Bleu, — Abu) I 


after a little algebra, using A/B = X. 
It is plain from (2) that the sign of Ad is the same as that of BAA — AAB. Now for 
the first step in the process, \ = X, , and hence from (5) it is obvious that 


AX < 0 if Ay < Onn / Der ’ 
(6) 
Ar > Of Ay > aee/Dex 


This leads to the criteria for convergence of the relaxation method: 

(i) in order that the process may converge to the lowest mode, the initial values z 
should be chosen so that the initial value of X, is less than a,,/b,, for all k; 

(ii) in order that the process may converge to the highest mode, the initial values z 
should be chosen so that the initial value of \, is larger than a,,/b,, for all k. 

It will be explained shortly how these choices of the initial values x can be made. 

If the choice (i) is made, the first step of the relaxation process, with \, = X, leads 
to a new value x with a smaller \. If \, is changed to this smaller \ and the process 
repeated, it is plain that the value of will decrease further. In practice it is not con- 
venient to change , after each step in the process, and if this change is not made we 
cannot use (5) to argue directly that a further step will decrease ) still further, for the 
d in (5) would be less than \, . The sign of the first term in the bracket in (5) will always 
be positive, but the sign of the term involving (A — \,) may be positive or negative. 
However, if \, is near to the lowest proper value the value of (A — A,) will be small, 
since from Rayleigh’s principle \ must always be greater than the lowest proper value. 
Even if \, is not near to the lowest proper value, (A — X,) will be small if the process 
is not carried on for too many steps with a fixed \, . In any case, no matter how many 
steps of changing an x, to reduce an (A, — A,B,) to zero are carried out with a fixed 
X, , A will not increase above this \, ; for (5) shows that as soon as \ approaches \, its 
value will begin to decrease. It is therefore plain that the process will lead to continually 
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decreasing values of the \ calculated after a series of steps. In practice the relaxation 
process is continued with a fixed value \, , so long as the residuals (A, — ,B,) as a 
whole can continue to be decreased. It should also be remarked that for any assumed 
A, , one solution is given by x, = 0 for all k; the computer must avoid approaching 
this solution by multiplying each x, by a factor from time to time so that the greatest 
of the x, is kept at some constant magnitude. 

Similar remarks, with obvious modifications, apply if the choice (ii) is made. \, 
then increases continually. 

Since the \ are bounded above by the highest and below by the lowest characteristic 
values of the equations, the \, must tend to a limit. This limit must be a proper value 
of the system, for all the equations (1) are then satisfied. The proper value in question 
could conceivably be one other than the highest or the lowest, but if this were to occur 
in actual computation it could do so only as a result of rare good fortune, and by choosing 
a higher (or lower) value of \, than the proper value thus found we could proceed to 
find the actual highest or lowest proper values. 

2. Practical applications. We shall now show how a suitable set of values x to serve 
as the first approximation in calculating the highest mode can be found. The method 
applies equally well, with an obvious change, to the lowest mode, but in most physical 
problems a good approximation to the lowest mode can be guessed intuitively. Let r 
and s be the values of m for which @,,,,/Dnm takes on its highest values. The highest value 
of X possible with all z,, zero except x, and z, is given by the larger root of 


ay, — Ab., ay. — Ab... 
= 0, 
ay. — Ab,, ay. — Ad,, 


a quadratic which is easily solved. The corresponding values of x, and x, are given by 
(a,, Po Aob,,) 2, + (ay, = hb,,)2, _ 0, 


where ), denotes the larger root of the quadratic. These values of z, and z, , with the 
other x’s zero, should be used as the starting point for the relaxation process. 

3. Estimation of frequencies intermediate between the highest and lowest. In the 
relaxation process described in the references below, it is necessary to find the lowest 
mode first, and then in finding the other modes to correct the values of the z’s found 
in the successive steps so that they continually satisfy 


> a,,.252, = 0, . (7) 


where x‘"? denotes the value of x, in the lowest mode. The fact that the highest mode 
can be found directly saves some time in this process. We shall now show that after finding 
the highest mode we can alter the problem to a new one so as to make the next highest 
mode of the old problem the highest mode of the new problem. When this is done, the 
steps for finding the other modes can be carried out without corrections to satisfy (7). 

Let \“? denote the i-th proper value, z,"? the corresponding mode values of x , 


then we have 


Lome ox” 2 ht". (8) 
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From this it follows by a well-known argument that if \“” ¥ »“” 
2s Gets 2G” = 2 b,25° xh” = @, (9) 


and even if two proper values are equal we can choose the modes corresponding to them 
so as to make (9) true. We shall also suppose the modes normalised so that 


( (4) 
Dd a,,.26?2” = 1. (10) 
rs 
Now consider the matrix A‘ = | a}! |, where 


al? wD at” >. «.at”. (11) 
A e 


We have 


Il 


Dd ala? = Y anzk? SY a,.26? 2h 
8 k es 


( 0 if ij, (12) 


» G2, if t=). 


Hence, from (8), 
( 0 if i #1, 
> (a. — af?)zi? = 4 
. a Fs ee if jt. 


After finding the highest mode and proper value, say \‘"’, we can form the matrix 


1 (n) 
ars — a;s — G,, . 


The highest mode of the original system of equations is annihilated by this matrix, but 
in the system of equations 


> az, =r > b,.2., (13) 


all modes save the highest remain proper modes with their former proper values: while 
the former highest mode now corresponds to the proper value zero. On finding the highest 
mode of the system (13) by the method described above, we get the second highest 
mode of the old system. This new mode can then be eliminated in its turn and the mode 
below it found. 

The procedure recommended, therefore, is to find the lowest and highest modes in 
the normal manner, and then to find the second highest, third highest modes, etc., in 
order, in the manner described. 

It is clear from (12) that the matrix |a,, — te as; | annihilates all the modes; 
and since these modes are independent vectors and every vector can be expressed as a 
linear sum in terms of them, it must annihilate every vector and therefore must be the 
zero matrix. We have therefore 


n 
(4) 
ie 2a. 


t=1 
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After (n — 1) modes have: been found, and the reduction of the matrix a,, carried 
out for each of them in the manner described, we must be left with the matrix 


(i) (4) 
ars — > Gre = Ay, ’ 


ima 


where 7 is the number of the mode left to be found. In this matrix the terms of any row 
are proportional to >>, a,,x‘"” so the remaining i-th mode can be found by solving the 


set of linear equations 
(i) 
bs A,X, —_ a, ° 
8 


The reductions can often be carried out with advantage in terms of the matrix b,, ; 
this is done by normalizing the modes so that instead of satisfying (10) they satisfy 


> ; b..2 i? _ 1, 
re 
and calculating aS! from the formula 
(@) (a) (i) (4) 
GS = A P bm bY pS breYe ) 
k e 


which is easily deduced from (8). This is particularly convenient in the most usual type 
of problem in which 6,, is a diagonal matrix—in mechanical problems, those in which 
the kinetic energy can be expressed as a sum of squares. If | b,, | = | m,é,, | the modes 
must be normalized by 


> my’ = 1 


and then 
(t) } (i) (4) 
a? = XO 'm.m,x? 2. 


After reductions corresponding to all but the 7-th mode have been made on the matrix 
a,, , the remaining matrix is a‘!’ from whose rows the 7-th mode can be found imme- 


diately. 


A NORM CRITERION FOR NON-OSCILLATORY DIFFERENTIAL EQUATIONS* 
By AUREL WINTNER (The Johns Hopkins University) 


Let f(é), x(t), A(@), --- denote real-valued, continuous functions on an unspecified 
half-line, 4; S t < @. If X(é) is positive on this half-line, put 


A* = r*() = AH) | (du)/2*(u), (1) 

t 
provided that the second factor on the right of (1) is a convergent integral. Under this 
proviso, a direct substitution of (1) shows that, if A¢¢) is a solution of the differential 


equation D,(A) = 0, where 


*Received Nov. 17, 1947. 
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DX) = DAD) = WO +IOAO, = (’ = d/dt) (2) 


then D,(A*) = 0, i.e., that A*(é) represents another (linearly independent) solution. 

Following A. Kneser’, let the differential equation D,(x) = 0 be called oscillatory 
or non-oscillatory according as each or none of its solutions x(#) ¥ 0 has zeros clustering 
at 1 = o. This alternative is complete, since, in view of Sturm’s separation theorem, 
either every or no solution z(t) # 0 of D,(x) = 0 has an infinity of zeros on the half- 
line i = t . The decision of the alternative (for a given coefficient function, f = f(t), 
of D,) is fundamental in certain questions of stability and related applications’. 

It seems to be of both theoretical and practical interest that the decision can always 
be based on a criterion similar to the “norm” conditions in the theory of linear func- 
tionals and operators (Lebesgue-Toeplitz). It is a criterion the applicability of which 
does not involve, in principle, the knowledge of a solution x(t) # 0 of D,(x) = 0, since 
it depends on the consideration of arbitrary functions. It can be formulated as follows: 

The differential equation D,(x) = 0 is of non-oscillatory type if and only if there exists 
some positive function, say X(t), corresponding to which the assignments (1), (2) define 
two continuous functions the product of which is absolutely integrable, 1.e., 


[ d* | D(A) | dt <@. (3) 


As an illustration of how to apply this criterion, choose the arbitrary function X(t) 
to be t. Then (1) and (2) reduce to A* = 1 and D,(A) = f(é)t, respectively, and so (3) 
will be satisfied if | f(t) | ¢ has a finite integral over the half-line. It follows that the 
absolute integrability of f(i)t (which, incidentally, is compatible with lim sup f(é)t = 
co and lim inf f(t) = — ©, where t— ©) is sufficient in order that the differential equa- 
tion D,(x) = 0 be of non-oscillatory type. 

Actually, this particular sufficient condition is contained in an asymptotic result of 
Bécher®. But this is not a necessary condition. In fact, other sufficient conditions result 
if the choice \ = ¢ is replaced by other choices of the arbitrary. function \(t). Such 
choices can be made relative to the coefficient function, f, of D,; , rather than in a way 
which, as in \ = #, is independent of f. 

Proof of the sufficiency. This part of the italicized criterion can be deduced from the 
following fact, which is a corollary of a general theorem*: If p = p(t) ~ 0 and gq = 
q(t) are continuous functions for large positive ¢, then the condition 


[10 1( J" 10 [7 au) at <2 (4) 
is sufficient in order that some solution y = y(t) of the differential equation 
(py +q=0  ('=d/dt) (5) 


1A. Kneser, Untersuchungen tiber die reellen Nulistellen der Integrale linearer Differentialgleichungen, 


Mathematische Annalen 42, 409-435 (1893), p. 411. 
2T. v. Karman and M. A. Biot, Mathematical methods in engineering, New York and London, 1940, 


Chapter VII and the references on p. 322. 
8M. Bocher, On regular singular points of linear differential equations of the second order whose coeffi- 
cients are not necessarily analytic, Transactions of the American Mathematical Society 1, 40-52 (1900), 


pp. 48-52. 
4A, Wintner, Asymptotic integrations of the adiabatic oscillator in its hyperbolic range, to appear in the 


Duke Mathematical Journal 15, (1948). 








PIES ne 


aes 











1948] AUREL WINTNER 185 


should tend to a finite limit, as ¢ ~ ©, and that this limit, y(@), be distinct from 0. 
It follows that, if \ = X(t) is any positive function possessing a continuous second 

derivative, then the case 
p =X’, q= dd" + fr (6) 


of (5) must have some solution y = y(t) which does not vanish from a certain ¢ onward, 
if condition (4) is satisfied by the functions (6). But it is clear from the definitions (1), 
(2) that the case (6) of (4) is identical with (3). Since X(t) is positive, it follows that (3) 
implies the non-oscillatory character of that differential equation for x = x(t) which 
results when y = 2/) is substituted into the case (6) of (5). 

The result of this substitution is seen to be the differential equation 


(x’n — xr’)’ + (A + fra = 0. 


Since the latter can be contracted into (x’’ + fz)A = 0, where A > 0, it is equivalent 
to x’ + fx = 0 and so, in view of (2), to D,(x) = 0. This completes the proof of the 
sufficiency of (3). 

Proof of the necessity. This part of the criterion is of theoretical interest only, and its 
verification is straightforward. As a matter of fact, \(¢) can now be chosen to be a solution 
x(t) of D;(x) = 0. 

In order to see this, suppose that the differential equation is. non-oscillatory. Then 
there exist a constant f, and a solution z(t) of D;fr) = 0 such that z(t) > 0 when t S 
Let ¢° be any value exceeding ¢, , restrict ¢ to the half-line ? < t < », and put 


< t 


A(t) = x(t) [ (du)/x*(u). (7) 


Then A(é) is positive, since x(t) is. Furthermore, it is easily verified from (2) and (7) 
that D,(A) = 0, since D,(x) = 0. Hence, in order to prove that condition (3) is satisfied 
by the function (7), all that remains to be ascertained is that the function (1) exists 
in the case (7), i.e., that 


i (du)/(u) < @ (8) 


holds by virtue of (7). But this can be ascertained by an elementary argument used 
by Hartman’. 

In fact, it is readily verified from (7) that the Wronskian, x\’ — Az’, of x(t) and 
\(2) is the constant 1. Hence, the derivative of x/) is identical with —1/’, and so 


x(t)/A(t) = const. — | (du)/d*(u). 


“to 


Since x(t) > 0 and X(t) > 0, it follows that 


0 < const. — | (du)/d*(u). 


This proves (8). 


5P, Hartman, On differential equations with non-oscillatory eigenfunctions, to appear soon. 
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THE EXTRUSION OF PLASTIC SHEET THROUGH 
FRICTIONLESS ROLLERS* 


By G. F. CARRIER (Brown University) 


1. Introduction: The Saint Venant-Mises theory of slow plane plastic flow has 
repeatedly been applied to problems concerning the deformation process which occurs 
when sheets are formed by passing the material between two fixed cylindrical surfaces 
with parallel axes. These problems include, of course, the problems of the rolling, ex- 
trusion, and drawing of such sheets. The analyses given are of two kinds. In one [1],’ 
[6], a simple one-dimensional theory is given; in the other [2], a more laborious scheme 
is used wherein the flow field is determined numerically by the method of characteristics 
(now very familiar to engineers because of its use in supersonic aerodynamic theory). 
It seems desirable to present information which either justifies the use of the simple 
theory, or modifies the simple theory so that its results become as accurate as those 
obtained in the numerical process mentioned above. When the thickness ¢ of the formed 
sheet and the diameter R of the cylindrical forming surfaces are of the same order of 
magnitude, it appears that the numerical scheme mentioned above is the most efficient 
procedure. However, when the sheet is thin (¢/R < 1) this procedure becomes very 
tedious. In this paper, we develop, from the fundamental equation of the Saint Venant- 
Mises theory, an approximation technique which leads directly to a justification of the 
one-dimensional theory for the cases where the cylindrical surfaces are frictionless and 
t/R <« 1. The less idealized case will be treated in a subsequent paper. 

2. Formulation of the problem: As is well known [3], the stress analysis in the case 
of problems of plane plastic strain is based on the yield condition 


(co, — o,)? + 4r:, = 4k’ (1) 


and the equations of equilibrium 


0a Oo 
— ——— = 0), 2 
Ox Oy (2) 
0a OT 
—< —— = 0). (3 
Oy Ox } 


The following substitution of variables simplifies the procedure. We set 


o,/k = 2w + sin 20 (4) 
a,/k = 2w — sin 26 (5) 
Tz,/k = — cos 20. (6) 


We note that this manner of expressing the stress components implies the satisfaction 
of Eq. (1). 


*Received Nov. 24, 1947. The results presented in this paper were obtained in the course of research 
conducted under Contract N7onr-358 sponsored jointly. by the Office of Naval Research and the Bureau 


of Ships. 
1Numbers in brackets refer to the bibliography at the end of this paper. 
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Omitting details which can be found in Reference [4], we observe that under fairly 
general conditions Eqs. (2) and (8) are equivalent to the system 


Ye + x, cot 6 = 0, (7) 
Ye — x, tan 6 = 0, (8) 


where 
E=ot 8, n=w- 6. 

It is evident that this system of equations is linear in the functions xz, y which are to 
be determined as functions of — and 7. It must be noted that this transformation can 
be applied only when the region in the x, y plane does not correspond to a line or point 
(degenerate region) in the £, 7 plane. This will occur in a limited part of our flow region 
[4]. 

The boundary conditions are related to the stress conditions at the cylindrical surface 
and at the inlet and exit sections. In the absence of friction the shear stress tr must vanish 
at the roll surfaces. On each of these surfaces then 


—cos [2(6 — y)] = 0, 


T(X, Yo) 


tan y = dy,/dz, 





exit section 














Fig. 1. 
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where y¥ = Yo(x) is the equation of the roll surface. If we consider that the material 
enters and leaves as a rigid body (i.e. if we assume the elastic deformations before en- 
trance and after exit to be negligible), we must demand that the material leave the roll 
(see Fig. 1) with uniform horizontal velocity. Since the velocities must be continuous, 
a uniform state must then exist in the plastic region adjacent to exit. In fact, it is known 
(see, for instance, Reference [2]) that region I (Fig. 1) is associated with single values 
of — and 7 or of w and @, namely: w = w,,,, = w, , 9 = 2/4. The former number is asso- 
ciated with the pull exerted externally on the exit section and the latter value follows 




















from the condition that r,, must vanish on y = 0. It is also known (see for instance 
Reference [5]) that since region II contains a straight characteristic (its boundary with 
region I) it is a region in which = const. = 7/4 + w, , but where 7 varies.” The map- 
ping of the z, y plane onto the é, » plane is shown qualitatively in Fig. 2. 
\7 
\ 4” oa 
\ Pg g= 7/4 
a Fe rd 
+ § 
“IN P 
% 
Ill Q 
Vx 
mk 
/ -* 
/¥=Yo(*) 6 
Fig. 2. 


in view of the degeneracy of the mapping of regions I and II and in view of the 
previously stated remark on such mappings, it is evident that any approximate solution 
of Eqs. (7) and (8) could not apply to regions I and II. We should also note at this 
point that a boundary value problem of this type cannot have a single analytic solution 
for the entire region. In fact, the derivatives £, , £, , 72, 7, Will, in general, be discon- 
tinuous at the boundaries separating the individual regions. Thus, zt zs not possible to 
write with rigor either 


P= Lasie,y, 1= LD bslery). 
r= Daw, y= Dba 2) 
or any other such development. 
However, we can write 
t= > agpi(E, 0) = 3 [g:(w)](@ — 2/4)" (9) 
y= > Byi(E, 0) = > [f(w)](@ — /4)' (10) 


‘ase where non-vanishing wall friction exists, but it 


?This remark would be slightly modified for the 
applies rigorously to our present problem. 
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where by proper choice of N, the f; , and the g; we can make these functions approach as 

closely as we wish* the functions x(w, 0) and y(w, 6). We observe that the actually existing 

discontinuous derivatives will not appear in these functions, but we also note that no 

differentiation of these functions is required in order to give the complete state of stress. 
The equation of the upper roll boundary can be written as 


y(x)/R =e+1-— (1 — 2°/R’)'”? =€+ 1 — cosy. 
On the roll surface (i.e. at y = Yo) 


7 = —cos 2(6 — y) = 0, 


6=n/4+7, (11) 
and on y = 0, 7,, vanishes; that is, on y.= 0, @ = 2/4. We now let the quantities z, y 
in Eqs. (9) and (10) be the dimensionless coordinates* +/R, y/R, so that on the roll 


boundary we have [from Eq. (11)] 


6 = n/4 + 7 
y=e+1-— cosy, 
x= siny, 
that is, 
N N ’ 
siny = >> gi(w)(0 — 2/4)’ = Do gilwy’, (12) 
N , 
e+1— cosy = > filwh’. (13) 


We now note that Eqs. (7) and (8) can be written 


(1 + tan y)y; + (1 — tan y)z; = 0, 


(1 — tan y)y, — (1 + tan y)z, = 0, 
or 
(l+y+7°/3 + --:)¥et (lL —7 t+ +°-)m = 8, (14) 
(l—-y—-<::)y—-Q+tv+--:)a, =9, (15) 
and when we use Eqs. (9) and (10), we obtain [noting that y(w) is odd in @ — 1/4] 
fi _ — 9 ’ (16) 
J2 = go + go'/2, (17) 
, (18) 


f= 5 (295 + gt’ — gh’’/2), 





3In the integral mean square sense. 
‘R is the radius of the cylinder. 


- 
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1 (248 = ie oe te") ( 
i, = 4 3 + + 3 ’ (19) 


When these are substituted in Eqs. (14), (15) we obtain two differential equations in 
the functions y(w), go(w). These equations can be written in the form 


ee ie, J rs he = (4 ra? , 1 7 
V2+et 1 Gate \\3% F 3H — 3% 


. - we ¢ De. B wis 
+ (206 — 459 — 239)’ + 6% + dsr (20) 





1 ee Mas. l wv = 
+ 2( + 99 yy’ Y Tt (245 — 39 + 39 + 390 (1 vy’ + \ 


mH | I 2 
— Bodo _ +i= _— oo (295 + w5)( 06 + + 59 \y 


go/e+ 1 2(go/2 + €) 
nos Lay)» — (2 : 
(dag + 67° + 3% 4295 + Jo) 9s — 390 (21) 
Pare Lag — gr 4 gt — Bah 4 
F 39 1: 3 hy + (Bos — 35" — 780% 12% — 75%) + 


We note that these equations are of an order which depends on the value of N [see 
Eq. (12)]..Thus, many boundary conditions are indicated. Recalling, however, the given 


information concerning regions I, II, we see that at 6 = 2/4 [i.e. y = OJ, go(w.) = «. 


Also at y = 2e (the upper corner of region II), we have w + 0 = w, + 1/4, 0 = 1/4 + 








2e hence w = w, — 2e. If we now define s = w, — w, the differential equations (20), 
(21) become 
a. = |] wee eee {($00 4 2 = Lily 
e+7/2~ (7/2 + 6 \\3% © 3% — 3% 
2 ] 1 iv ] v 5 
Se. Fst ese sae =. 5 9 
5 (246 459° 2396 + 62° + sip (22) 
1 ” 2° 2 , 4 + ] wr 1 iv ’ 
+ 2( + 99 0 by ¥ + (246 = 390 + 390 4 + 37 yy +t: ae 
a be r ; si) 
Get 1 ~*~ 22+ | OMT WH + of 
— (95 + Laer + dyer? — Hogs + oi( 20a — 45" (23) 
39 RE 39° Jo Jo Jo 39° e 


Love tmp fh 2 Bn Sa Rely wf 
+ 308” + 50s hat + (Toot — gg" — sani” + aoe” — age + of 


where the primes now indicate differentiation with respect to s. 
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Rigorously, we have more boundary conditions along the line § = const. which 
represents region II of the physical plane. If these were applied, the problem would be 
over-determined and, in fact, we could at best obtain a solution valid for region III 
only. If we relax these conditions and define g = In(1 + y’/26), ¥ = In(1 + gi/26), 
Eqs. (22) and (23) can be put into the form 


= o—_ ¥ - 1 —— 4 4? 2 7 1 ) 3 
re | (7/2 + ©) {(§o6 F 39 — 390 Jy 
2 1 1 iv 1 :) E 
oie i —_—s Q,/079 rca —~. 5 
(306 459° — 2390" + G90 + Gogo} (24) 
l yz ie 2 , 4 ”r ] tr 1 iv , 
2 Aly, + 990 )y y + (206 = 3% + 3% + 39° Jay + -: : ds, 


go’ 


_ | ; gi’ \ » 
y=s+t+ | Dg2/2 + ©) {206 + oid + 2 by 


I ? l vt? 3 1 , 2 , 4 r 
Jo + 59 hy + 4296 + 0.)(296 _ 39 (25) 


1, 
. (doc t 5 3 


+ 08” ~ Lay - (Svs — z508" — Sah" + a0 — dy" + a ds 

The equations are of the conventional non-linear Voltera type and may be treated 
by successive approximations. Without going into details, it may be stated that the 
first approximation and the subsequent ones are in excellent agreement for s > 2e, 
whereas in the neighborhood of s = 0, the behavior is fairly erratic. However, y(s) is 
of interest only for s > 2¢ and go(s) only for s > e. Hence, the solution 


y(s)~s, o(s)~s (26) 


provides an excellent estimate of the state of stress existing in the flow field. We can 
write in fact, [this is equivalent to Eq. (26)] 
w, — w = In(l + 7’/26) (27) 
along the boundary, and 
w® —w = In(1+ 9 Qe) = In(1 + 2°/2¢e) (28) 
along y = 0. These formulae are good for s > 2e and are essentially Sachs’ formulae [1]. 
The corresponding results for the case where the wall friction does not vanish require 
certain numerical work and will be presented in a later paper. Detailed interpretative 
remarks concerning the foregoing result do not seem to be in order here since they would 
necessarily coincide with the findings in [1] regarding this specific problem. 
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ADDITIONAL CORRECTIONS TO OUR PAPER 
THE CYLINDRICAL ANTENNA: CURRENT AND IMPEDANCE 
QUARTERLY OF APPLIED MATHEMATICS 3, 302-335 (1946) and 4, 199-200 (1946) 


By RONOLD KING anp DAVID MIDDLETON (Harvard University) 


page 305, Change the number of Eq. (13a) to (13); delete “‘so that” following Eq. (13a); 
delete Eq. (13b). 
page 306, Change the number of Eq. (14a) to (14); delete “where” following Eq. (14a); 
delete Eq. (14b). 
Eq. (16)—add superscript —1 on R,, in the integrand. 
page 317, Fig. 10—The value | ¥.(h — \/4) | should be at 16.6 instead of 17.4 with 
appropriate changes in the several curves. 
page 321, Figs. 12 and 13—All the curves are somewhat in error for Bh < 2/2. The 
correct values are obtained from (74), using the corrected values for ¥ obtained 
from Fig. lla on page 200 of volume 4. 
page 327, Table II—First line: Insert — between 2/2 and Bh,,, . 
Second line: Replace 800 by 820. 
Fourth line: Replace 67 by 73. 


page 328, Eqs. (14a), (14b)—Insert — after =. 


page 330, Eq. (23)—-Change sign of lower limits on all three integrals by inserting — 
sign. This change is in addition to corrections on page 200 of volume 
4, 
Eq. (27)—Change first — sign to +; change last + sign to —. 


page 335, Eq. (45)—Change all upper limits in four integrals to u, . 
Change all lower limits in four integrals to —1u, . 
Eq. (46)—Last integral only: Change upper limit to u, , lower limit to —1, . 
Eq. (47)—Delete superscript bars in second integral of the first member of the 
equation. Change — to + before this second integral. 


*These articles have been translated for the David W. Taylor Model Basin, United States Navy, 
by the Applied Mathematics Division, Brown University, Providence, R. I. 














